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Little is known of the internal structure or composition 
of Mars or the present rates and characteristics of tectonic 
activity and meteoroid impacts. The scientific tool best suited 
to address these issues is seismology. While a simple seismic 
experiment was included on the two Viking landers, one 
of the instruments did not operate and the experiment was 
severely limited by the location of the sensor high on the 
lander, the low sensitivity of the instrument, and the limited 
data rate allocated to seismic measurements. Thus, even 
after the Viking mission, our knowledge of the seismic char- 
acteristics of Mars remains poor. 

The Mars Global Network Mission, now in the early stages 
of mission planning, offers an opportunity to deploy a global 
seismic network on Mars for the purposes of determining 
the internal structure and constitution of the planet and 
the characteristics of marsquakes and meteoroid impacts. 
At the request of the Mars Science Working Group, a work- 
shop was held May 7-9, 1990, in Morro Bay, California, 
to define the scientific rationale and technical requirements 
for a global seismic network on Mars. This report is a sum- 
mary of the findings and recommendations of that workshop. 

The principal scientific objectives of a martian global 
seismic network are (1) to determine the nature and struc- 
ture of the martian crust; (2) to determine the structure 
and state of the martian mantle; (3) to determine the radius 
and state of the martian core; (4) to determine the locations 
and mechanisms of marsquakes; (5) to determine the impact 
flux and characteristics of meteoroids in Mars-crossing orbits; 
and (6) to provide supporting information for meteorology. 

From the experience of the Apollo seismic network, it 
is clear that when we design and conduct a seismic exper- 
iment on a new planet, (1) we must not assume that seismic 
signal characteristics are similar to those we see on Earth, 
(2) we must design an instrument with the highest possible 
bandwidth, dynamic range, and sensitivity, and (3) we must 
be alert to the possibility of unexpected seismic sources. 
The Viking seismic experiment, even though it was highly 
curtailed by mission priorities and by the inability to uncage 
the seismometer on the Viking 1 lander, provided valuable 
information for the planning of future missions to Mars. 
First, it established that the seismic background noise on 
Mars due to winds and atmospheric pressure fluctuations 
is very low. Seismometers more sensitive than the Viking 
instrument by a factor of at least 10 4 can operate on the 
planet without being affected by typical martian winds. 
Second, a much broader frequency response is required than 
was used for the Viking seismometer, opening up the possi- 
bility of detecting normal modes, surface waves, tidal loading, 
and Chandler wobble. Finally, a significant network of instru- 
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ments will be necessary to study seismicity and to deter- 
mine internal structure. 

The principal natural seismic events on Mars are mars- 
quakes and meteoroid impacts. While the rates of occur- 
rence are not known for either type of event, theoretical 
considerations and experience from the Moon provide a 
basis for their estimation. Thermal stress associated with 
the global cooling of the martian lithosphere is expected 
to give rise to marsquakes at a rate of more than 10 events 
with seismic moment in excess of 10 23 dyn cm per Earth 
year and more than 2 events with seismic moment in excess 
of 10 24 dyn cm per Earth year; events in this size range 
are expected to produce records observable throughout a 
global network. Additional possible contributors to litho- 
spheric strain rate (and thus to expected seismicity) on a 
regional to global scale include enhanced cooling beneath 
major volcanic provinces (e.g., Tharsis), seasonal variations 
in polar cap loading, and membrane stresses induced by 
tides and by long-term changes in planetary obliquity. 
Potential sources of marsquake activity on a more local 
scale include magma motion, landslides, variations in solar 
insolation, and freeze-thaw cycles. We conclude that internal 
seismic activity will occur on Mars at a rate more than 
adequate to address the principal scientific objectives. 

Available data indicate that we can expect impacts of 
both cometary meteoroids and asteroidal fragments at rates, 
per unit area, on Mars similar to those detected on the 
Moon by the Apollo network, with the impact rates of 
asteroidal objects possibly somewhat higher on Mars if there 
is an abundance of meteoroids with Mars-crossing orbits. 
Most meteoroids of cometary origin are likely to be effectively 
consumed in the atmosphere because of their low density, 
friability, and high encounter velocity, while those of 
asteroidal origin in the mass range that yielded detectable 
seismic signals on the Moon are likely to impact the surface. 
Overall, given the greater surface area of Mars than the 
Moon, but also the likely higher seismic attenuation in the 
interior, the effect of the martian atmosphere, and the 
probable intervals of significant wind-generated ground 
noise, the teleseismically detectable seismicity rate from 
meteoroid impacts of asteroidal origin is likely to be some- 
what less on Mars than on the Moon, and that due to 
cometary impacts is likely to be substantially less. 

For each of the identified scientific objectives, certain 
classes of seismic observations are likely to provide the most 
straightforward means to achieve the objective. The need 
to make these key observations in turn affects station siting 
plans. With a global network, crustal structure is most readily 
addressed by measurement of the dispersion of 10-100-s 
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surface waves; such a measurement can be made with a 
small array of three stations, even for sources with poorly 
known locations and origin times. Mantle structure is most 
readily addressed by the inversion of the travel times of 
body waves vs. epicentral distance; in particular, major dis- 
continuities at depth (e.g., the olivine-/? phase boundary) 
should be detectable through pronounced triplications in 
the travel-time curve. The radius and state of the core 
will be inferable from the measurement of the travel times 
and amplitudes of core-reflected body waves at short 
distances (0-20°) and from the location of the core shadow 
zone and the arrival times of core phases at large distances 
(100° and beyond). It is thus important that some seismic 
stations be located at considerable distance from areas of 
potential seismic activity. The state and structure of the 
core strongly affect the periods of the gravest free oscillations 
of the planet, but these modes (with periods as great as 
2000 s) are not excited except by very large seismic sources. 
The location and characterization of natural seismic sources 
is aided by having seismic stations near the source; such 
a consideration does not affect station siting for effectively 
random sources (e.g., meteoroid impacts), but does call for 
the placement of several stations near areas of more likely 
tectonic seismicity (e.g., Tharsis). The objective of providing 
information in support of meteorological experiments (e.g., 
on surface boundary layer parameters) suggests that seismic 
and meteorological stations be cosited where respective siting 
criteria permit. 

There are technical trade-offs for seismic experiment 
performance between emplacement of the seismic sensor 
beneath the martian surface (i.e., using a penetrator) and 
placing the sensor package on the surface from a lander. 
The subsurface package has superior coupling to ground 
motion, is far less susceptible to surface temperature 
variations and wind stress, and should yield a higher signal- 
to-noise ratio. The current lack of availability of a long- 
term power source sufficiently rugged to withstand the 
decelerations of penetrator impacts, however, apparently 
precludes consideration at present of subsurface emplace- 
ment for stations in a long-term network. Surface stations 
should be adequate for the recording of seismic signals at 
periods shorter than 25-100 s, which covers the desired 
bandwidth except for the longest period surface waves, 


normal modes, and tidal and solid-body excitations (e.g., 
Chandler wobble). It will be essential, however, to site the 
seismic sensors as far as possible from the lander, which 
is likely to be the principal local source of seismic noise, 
particularly if the lander presents a large cross section to 
near-surface winds. The sensors should also be installed in 
streamlined enclosures. 

Other desirable specifications for a seismic station on Mars 
are as follows: number of data channels - 3 components 
of ground motion; A/D sample rate - 1000 sps; A/D -24 
bit; stored data rate - 50 sps; total data rate - 100 Mbits/day; 
bit error rate -<10"*; bandwidth - <0.04-20 Hz; sensitiv- 
ity - 10~'° g; mass - 1.5 kg (sensors only); volume - 0.01 m 3 ; 
power - 2 W; onboard RAM - 4 Mbytes; calibration - once 
per day. The total data rate presumes a simple data 
compression scheme that reduces data volume by a factor 
of 3. More powerful schemes can increase this compression 
factor somewhat, but triggering algorithms are considered 
scientifically undesirable in the absence of information on 
the types of signals likely to characterize martian seismic 
events. 

A station siting plan must strike a compromise among 
the various scientific objectives. There must be a mix of 
closely spaced stations to detect and locate nearby seismic 
events and to serve as arrays for measurements of phase 
velocity and propagation direction, and more globally 
dispersed stations to provide planetary coverage and to 
ensure the recording of core phases and other signals 
diagnostic of deep structure. The recommended plan consists 
of nested triads of stations. Small triads of stations should 
be located approximately 100 km apart as local arrays. Each 
of these small triads should consist of 3 three-component 
short -period instruments arranged in a triangle around a 
broadband observatory-class three-component sensor. Sets 
of these small triads should be emplaced in larger triangular 
patterns, approximately 3500 km on a side, one each in 
the eastern and western hemispheres of Mars. Tentative 
locations for large triangle vertices consistent with the 
various scientific objectives include (l)west of Ascraeus 
Mons, north of eastern Valles Marineris, and in the southern 
highlands and (2) in the Elysium province, in the Isidis Basin, 
and in the northern lowlands. 
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Introduction 


Little is known of the internal structure or composition 
of Mars. Also poorly known are the present rates and char- 
acteristics of tectonic activity and meteoroid impacts. The 
scientific tool best suited to address these issues is seismology. 
While a simple seismic sensor was included on the two Viking 
landers, one of the instruments did not operate and the 
experiment was severely limited by the location of the sensor 
high on the lander, the low sensitivity of the instrument, 
and the limited data rate allocated to seismic measurements. 
Thus, even after the Viking mission, our knowledge of the 
seismic characteristics of Mars remains poor. 

The Mars Global Network Mission, now in the early stages 
of mission planning, offers an opportunity to deploy a global 
seismic network on Mars for the purposes of determining 
the internal structure and constitution of the planet and 
the characteristics of marsquakes and meteoroid impacts. 
At the request of the Mars Science Working Group 
(MSWG), a workshop was convened to define the scientific 
rationale and technical requirements for a global seismic 
network on Mats. The workshop was held May 7-9, 1990, 
in Morro Bay, California, and was attended by 10 scientists 
with interest and expertise in seismology and martian geo- 
physics. A list of participants and the agenda for the 
workshop are included at the end of this technical report. 
This report is a summary of the findings and recommen- 
dations of that workshop. 

It is important at the outset to recognize that obser- 
vational seismology is quite different from meteorology, geo- 
chemistry, or most other planetary surface sciences. The 
nature of seismic signals and noise on a newly visited planet 
are unknown, and it is impossible to determine in advance 


all the optimum parameters of a global seismic experiment. 
In this respect, a seismic network is much like an imaging 
experiment. While we don’t usually think of a camera system 
as an experiment in this sense, the appropriate bandwidth, 
wavelengths, and dynamic range are also not known in 
advance. This lack of knowledge is generally accommodated 
by the enormous data rates and sophisticated processing 
that we have come to accept for imaging experiments. The 
lack of prior knowledge is a handicap only if weight, power, 
and data-rate constraints make it necessary to select narrowly 
defined experimental parameters a priori. Such a preselection 
is not done for imaging, and it should not and need not 
be done for seismology. 

Following a brief overview of the mission concepts for 
a Mars Global Network Mission as of the time of the work- 
shop, we present the principal scientific objectives to be 
achieved by a Mars seismic network. We review the lessons 
for extraterrestrial seismology gained from experience to date 
on the Moon and on Mars. An important unknown on 
Mars is the expected rate of seismicity, but theoretical 
expectations and extrapolation from lunar experience both 
support the view that seismicity rates, wave propagation 
characteristics, and signal-to-noise ratios are favorable to 
the collection of a scientifically rich dataset during the 
multiyear operation of a global seismic experiment. We dis- 
cuss how particular types of seismic waves will provide the 
most useful information to address each of the scientific 
objectives, and this discussion provides the basis for a strategy 
for station siting. Finally, we define the necessary technical 
requirements for the seismic stations. 
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The Mars Global Network Mission, as conceived at the 
time of the workshop, is to consist of about 20 landed stations 
globally distributed over the surface of Mars. These stations 
may be launched during a single Mars opportunity and 
deployed sequentially from orbit, or launched a few at a 
time over several opportunities (at approximately two-year 
intervals), with deployment occurring directly from hyper- 
bolic approach. There may be restrictions on allowable land- 
ing sites on the basis of altitude (altitudes ranging from 
-2 to 6 km are currently considered acceptable) and latitude 
(solar power is not thought to be viable poleward of about 
45°). An orbiting platform may be available for commun- 
ication support. 

The landers are envisioned as relatively simple systems, 
with a landed mass of about 75 kg and an electrical power 
system capable of supplying 5-15 W continuously. Decel- 
eration techniques will be used to minimize impact loads. 
Complexity of both mechanisms and operations will be kept 
to a minimum, and the stations should be capable of 
conducting operations on the surface for up to 10 years. 

The objectives of this mission have been chosen to take 
advantage of the unique opportunities presented by a global 
network of stations: the ability to make simultaneous 
measurements of a given phenomenon at widely spaced 
locations and the capability for sampling a large number 
of different geological settings. The primary objectives are 
scientific, and include global seismology, meteorology, and 
geochemical investigations. Other key objectives relate to 
future plans for the exploration of Mars, such as charac- 
terizing the upper atmosphere for optimizing aerocapture 
techniques and identifying potential sources of accessible 
water and other resources. 

In order to formulate a useful set of requirements and 
goals for a seismic network experiment on Mars, it is helpful 
to have a framework against which to judge their feasibility. 


Mission planners and engineers have developed several 
scenarios for a Mars network mission that attempt to address 
the technical constraints and science goals as they are 
currently understood. These scenarios have incorporated 
various assumptions about a seismic experiment, and result 
in a number of limitations as to how such an experiment 
could be conducted. Table 1 contains a list of parameters 
we consider relevant to seismology, taken from recent JPL 
and Ames Research Center network mission studies. While 
realizing that these values are subject to rapid change and 
may well be obsolete before this report is published, we 
include them here as our point of reference as to what 
is currently considered technically acceptable. In many areas 
we will argue for substantially more capability than is 
reflected in Table 1. 


TABLE 1. Initial planning constraints on a Mars 
global seismic experiment. 


Number of Stations 
Location Accuracy 

10-20 

Control 

50-200 km 

Knowledge 

<1 km 

Distribution 

Global 

Lifetime 

2-10 years 

Relative Timing Accuracy 
Telemetry Rate (per station) 

20 msec 

Via Orbiting Relay 

5 Mbit/day 

Direct to Earth 

<300 kbit/day 

Command Capability 

Little or none 

Emplacement 

On surface, <1 m from lander 

Landing Loads 

40-100 g 

Mass 

1.5 kg 

Power 

2 W 
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Scientific Objectives for a Global Seismic Network on Mars 


The principal scientific objectives of a martian global 
seismic network are to (1) determine the nature and struc- 
ture of the martian crust; (2) determine the structure and 
state of the martian mantle; (3) determine the radius and 
state of the martian core; (4) determine the locations and 
mechanisms of marsquakes; (5) determine the impact flux 
and characteristics of meteoroids in Mars-crossing orbits; 
and (6) provide supporting information for meteorology. 

DETERMINATION OF THE NATURE AND 
STRUCTURE OF THE MARTIAN CRUST 

Mars has a distinct low-density crust of variable thickness, 
as indicated by the partial to complete isostatic compensation 
of surface topography ( Phillips et at, 1973; Phillips and 
Saunders, 1975). The mean thickness of the crust, however, 
is poorly constrained. A minimum value for the average 
crustal thickness of 28 ± 4 km (depending on choice of 
crust-mantle density difference) was obtained by Bills and 
Ferrari (1978) by fitting a model crust of uniform density 
and variable thickness overlying a uniform mantle to 
topography and gravity expressed in spherical harmonics to 
degree and order 10; for this minimum mean thickness the 
crust is of zero thickness beneath the Hellas Basin. For 
a crustal thickness of 15 km at the site of the Viking 2 
lander, a result inferred from the tentative identification 
of narrow-angle Moho-reflected phases in the single recorded 
seismogram of possible tectonic origin (Anderson et al, 1977), 
Bills and Ferrari (1978) obtained a mean crustal thickness 
of 37 ± 3 km, a maximum thickness of 69 ± 8 km (beneath 
Tharsis), and a minimum thickness of 9 ± 1 km (beneath 
Hellas). In a contrasting work, Sjogren and Wimberly (1981) 
used topography and gravity data to infer that an isostatically 
compensated Hellas Basin has an Airy compensation depth 
of 130 ± 30 km. From the models of Bills and Ferrari (1978), 
this value would correspond to a globally averaged crustal 
thickness of about 150 km. 

The mean thickness of the martian crust bears strongly 
on the history of differentiation of the planet, in that the 
crustal volume consists of some combination of material 
remaining from postaccretional differentiation of crust from 
mantle plus later intrusive and extrusive magmas generated 
by partial melting of the martian mantle. A crust 30 to 
150 km in thickness constitutes 3% to 13% of the planetary 
volume, a substantial range of uncertainty. Seismology offers 
the only direct tool to measure crustal thickness, either 
by the travel times of body wave phases at near and regional 
distances or by the dispersive characteristics of surface waves. 

There are expected to be pronounced lateral variations 
in crustal thickness that bear strongly on the deep struc- 
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ture and mode of formation of major geological features 
on Mars. For instance, the global crustal dichotomy, the 
approximately hemispherical division of the martian sur- 
face between the topographically lower and stratigraphically 
younger northern plains and the heavily cratered southern 
uplands, is thought to be compensated by an Airy isostatic 
mechanism on the basis of analysis of gravity anomalies across 
the dichotomy boundary (Janie, 1983), but this hypothesis 
needs to be confirmed by seismic measurement of crustal 
thickness in each of the northern and southern hemispheres. 
The deep structure of the Tharsis province, which has been 
a focus for martian tectonic and volcanic activity over much 
of the history of the planet, has been a topic of considerable 
controversy. Structural models vary from a mix of Airy and 
Pratt isostasy (Sleep and Phillips, 1979, 1985; Banerdt et al., 
1982) to large-scale flexural support of the excess volcanic 
and intrusive material by the finite elastic strength of the 
martian lithosphere (Willemann and Turcot te, 1982; Banerdt 
et al., 1982); the distinct crustal structures predicted by 
these models can be readily distinguished with regional 
seismic data from the Tharsis area. In general, once the 
crustal structure is established seismologically in a few areas, 
extrapolation to the rest of Mars by means of gravity and 
topography data — both those now available and those 
expected to be provided by the Mars Observer mission — 
can be readily accomplished (e.g., Thurber and Solomon, 
1978). 

Significant reservoirs of H 2 0 and other volatiles are likely 
to be present within the martian crust, as permafrost and 
saturated layers of regolith and bedrock (e.g., Carr, 1987). 
The presence of ice and water at depth will significantly 
affect the seismic wave velocities and attenuation. While 
active seismic experiments offer the most promising approach 
to the exploration of buried volatile reservoirs, seismic phases 
recorded by passive stations at near or regional distances 
may nonetheless contain signatures of a layered crustal 
velocity structure indicative of ice-rich or water -rich layers. 

DETERMINATION OF THE STATE AND 
STRUCTURE OF THE MARTIAN MANTLE 

The only direct constraints on the internal structure of 
Mars are its mean density and its moment of inertia. The 
latter quantity is not precisely known, and, in the absence 
of a direct determination of the martian precessional 
constant, must rather be derived from the second degree 
zonal coefficient J 2 in the spherical harmonic expansion of 
the gravitational potential by means of one or more 
additional assumptions. The appropriate assumptions are 
matters of current debate (Bills, 1989a, b; Kaula et al, 1989), 

PRECEDING PAGE BLANK NOT FILMED 



6 Global Seismic Network on Mars 


and it may be concluded only that the dimensionless moment 
of inertia I/MR 2 most likely lies between 0.350 and 0.365, 
but lower values cannot be completely excluded. In addition, 
the interior constitution of Mars is constrained by geo- 
chemical and cosmochemical considerations relating the bulk 
composition to that of various classes of meteorites (e.g., 
Anderson, 1972), in particular by the hypothesis that the 
SNC meteorites were derived from Mars (e.g., McSween, 
1984; Bogard et al, 1984; Becker and Pepin, 1984). 

An early analysis of the plausible ranges in the density 
of mantle and core material and in the radius of the core 
was made by Anderson (1972), who concluded the mantle 
of Mars is richer in FeO than the Earth’s mantle. This 
analysis was updated by Goettel (1981), under the then 
commonly held view that the dimensionless moment of 
inertia of Mars was I/MR 2 = 0.365 ( Reasenberg , 1977; Kaula, 
1979). The zero-pressure density of the mantle under this 
assumption is between 3400 and 3470 kg/m 3 , confirming 
a magnesium number [molar Mg/(Mg+Fe)] somewhat lower 
than petrological estimates for the Earth’s upper mantle 
(Goettel, 1981). A comparatively FeO-rich composition has 
also been inferred for the source region of the magmas from 
which the SNC meteorites were derived (e.g., McSween, 
1985), a result that has been used as one of the arguments 
favoring Mars as the SNC parent body (Wood and Ashwal, 
1981). If the dimensionless moment of inertia I/MR 2 is as 
low as 0.345, as Bills (1989a) estimated, however, then the 
zero-pressure mantle density must be very low (3200-3220 
kg/m 3 ) (Bills, 1990), implying a high magnesium number 
and a difficulty with the identification of Mars as the source 
of SNC meteorites. While this controversy would be lessened 
by a precise measurement of the martian precessional con- 
stant, definitive determination of the magnesium number 
of the martian upper mantle must come from measurement 
of the mantle seismic velocity structure, either by an in- 
version of body wave travel times or by an inversion of 
normal mode eigenfrequencies. 

Determination of the radial structure of the martian 
mantle, of course, provides significantly more information 
than the mantle magnesium number. Prominent discon- 
tinuities, or rapid variations in seismic velocity with depth, 
can be signatures of solid-solid phase changes or of chem- 
ical layering resulting from whole-planet differentiation. For 
instance, the olivine-/? phase transition should occur at 
a depth greater than 1000 km (BVSP, 1981), with the pre- 
cise depth a function of mantle temperatures and magnesium 
number. Within layers of constant composition and phase, 
the radial gradients of compressional (P) and shear (S) wave 
velocities provide measures of the mean thermal gradient 
and thus constrain the efficiency of convection as a mantle 
heat transport mechanism. The nature of any asthenosphere 
on Mars — important for models of mantle dynamics, thermal 
structure, interior volatile budgets, and magmatism — will 


be revealed by zones of low seismic velocity and high seismic 
attenuation. 

Once the average radial structure of the martian mantle 
is reasonably well established, it may be possible to determine 
major lateral variations from that average structure. Such 
a possibility will depend on there being a fairly good global 
distribution of energetic seismic sources. Lateral variations 
in the seismic velocity structure of the lithospheric mantle 
could result from melting-induced compositional variations, 
perhaps most notably beneath major volcanic provinces such 
as Tharsis (Finnerty et al, 1988; Phillips et al., 1990). Lateral 
variations in the deeper mantle are likely to be the signature 
of patterns of mantle upwelling and downwelling convective 
flow (e.g., Schubert et al, 1990). 

DETERMINATION OF THE RADIUS AND 
STATE OF THE MARTIAN CORE 

While the moment of inertia of Mars indicates that the 
planet has a distinct high-density core, neither the size of 
the core nor its density are well constrained. For I/MR 2 = 
0.365, the core could vary in radius from 1300 to 1900 
km, corresponding to zero-pressure density varying from 
8900 (pure Fe) to 5800 (FeS) kg/m 3 and core mass varying 
from 13% to 26% of planetary mass (Goettel, 1981). For 
I/MR 2 less than this figure, the core radius is larger (and 
the mantle density less) for a given value of core density 
(e.g., Bills, 1990). From siderophile and chalcophile abun- 
dances in SNC meteorites, haul et al. (1986) and Treiman 
et al. (1987) have inferred that the SNC parent body has 
a core constituting 20-35% of the parent body mass and 
containing 12-14 wt% sulfur. The sulfur content of the 
core strongly influences the timing and extent of freez- 
ing of the core from an initially totally molten state; inner 
core freezing is a major potential source of energy for driv- 
ing a martian core dynamo (Stevenson et al, 1983; Schubert 
and Spohn, 1990). In particular, a totally molten core at 
present constrains the sulfur content (or that of other light 
elements in the core) to be at least 15% (Schubert and Spohn, 
1990). 

Seismic observations, including travel times and ampli- 
tudes of body waves and eigenfrequencies of longer period 
normal modes, are capable of determining the radius of the 
martian core, the radius of any solid inner core, and the 
average radial seismic velocity structure of the core. The 
size and average seismic velocity of the core will constrain 
core composition and thus the bulk composition of the 
planet, in particular the ratio of metal plus sulfide to silicate. 
Tie relative fractions of liquid and solid core will strongly 
constrain the fraction of light elements in the core, the 
thermal history of the core, and the availability of inner 
core freezing as an energy source for generation of a planetary 
magnetic field. 
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DETERMINATION OF THE LOCATIONS 
AND MECHANISMS OF MARSQUAKES 

The distribution and character of tectonic sources of 
seismic energy in a planet provide measures of the internal 
stress field and of the processes that have generated those 
stresses. Mars displays an abundance of tectonic features 
that must have formed by seismogenic faulting, and models 
for the present distribution of stresses within the martian 
lithosphere (e.g., Banerdt et aL, 1982; Sleep and Phillips, 1985) 
suggest that marsquakes should be generated wherever 
ongoing lithospheric processes (e.g., differential cooling, tidal 
or seasonal loading, fluid migration) induce further stresses 
that add constructively to the long-term stress field. 

Experience on the Earth and Moon suggests that a suitably 
configured seismic network on Mars operational over a dur- 
ation of several years will permit the determination of the 
principal regions of marsquake activity and should yield at 
least approximate source mechanisms indicative of the style 
of faulting. Such data will provide direct tests of current 
models for lithospheric stress on Mars. It can be anticipated 
that the most likely areas of seismicity are those that display 
geological evidence for relatively recent volcanic or tectonic 
activity, such as the Tharsis Montes, Olympus Mons, or 
Valles Marineris (e.g., Lucchitta, 1987; Tanaka et aL, 1988). 
Some densification of the seismic network in the most 
probable areas of seismic activity is warranted to ensure 
adequate resolution of marsquake hypocenters and source 
mechanisms. 

DETERMINATION OF THE IMPACT FLUX 
AND CHARACTERISTICS OF 
METEOROIDS IN MARS' 
CROSSING ORBITS 

Before the era of space missions, practically all the infor- 
mation about meteoroids, the small interplanetary objects 
in the inner solar system, came from observations of meteors 
and recovered meteorites. Various lunar and planetary 
missions have expanded this database to include crater 
statistics on most of the planets and satellites with solid 
surfaces, giving us information on the history and size 
distribution of impacts of relatively large meteoroids. 

The Apollo lunar seismic experiment, in particular, yielded 
a unique dataset that allowed determination of the flux 
and orbital characteristics of several types of meteoroids 
crossing the orbit of the Earth-Moon system (e.g., Oberst 
and Nakamura, 1991). Most prominent among these objects 


are those originating from the break-up of long-period comets 
and asteroids. Many of those of cometary origin are associated 
with meteor showers observed on Earth. In contrast, a 
majority of those considered to be asteroidal fragments are 
observed as sporadic impacts. However, some of them form 
streams and swarms, indicating relatively recent break-up 
from their parent bodies, and some of them even appear 
to be related to known meteorites. 

The true nature of these impacting objects is not clearly 
understood, and there are yet many open questions about 
their origin. Our present interpretation is limited by the 
fact that all our observations to date concerning the current 
population of these objects are restricted to those in Earth- 
crossing orbits. Similar observations on Mars will expand 
our knowledge to include those in Mars-crossing orbits. This 
will greatly advance our understanding of the nature and 
orbital properties of these objects, and thus will help us 
decipher their role in the recent evolution of the solar system. 
In particular, a determination of the relative impact fluxes 
on Mars and the Moon will provide a critical test of 
competing crater chronologies used for the dating of 
relatively young surfaces on Mars (e.g., Hartmann, 1977; 
Neukum and Hiller, 1981; BVSP, 1981). 

SUPPORTING INFORMATION 
FOR METEOROLOGY 

Seismic noise on Earth is generally of meteorological and 
oceanographic origin. Winds and waves couple into the 
ground and result in propagating seismic disturbances. On 
Mars this type of “noise” is actually a useful meteorological 
“signal” that is relevant to such boundary-layer parameters 
as surface friction velocity and surface roughness. Seismic 
instruments yield information about high-frequency wind 
parameters that are generally not measured by meteorological 
observations. 

Ideally seismic and meteorological observatories on Mars 
should be combined. The electronics and processing require- 
ments are similar and can be shared. Having the seismological 
and meteorological data corecorded helps in the analysis 
of both. The expertise of seismologists in time series 
processing can be used to advantage in designing the 
meteorological electronics, and the expertise of meteorol- 
ogists can be used to advantage in helping to interpret the 
background “noise.” Thermoelastic stresses due to prop- 
agating temperature gradients (sunrise, sunset, clouds, 
eclipses) and propagating weather fronts may also provide 
useful seismic sources for probing the shallow interior. 


Lessons from Apollo and Viking 
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There have been two previous experiments in extrater- 
restrial seismology, the Apollo seismic network that operated 
on the lunar surface for the 8-year period 1969-1977 and 
the Viking seismology experiment that operated on the 
martian surface for a 19-month period. Both of these 
experiments yielded valuable lessons for the design and 
operation of future seismic networks on other planets. 

THE APOLLO SEISMIC NETWORK 

The Apollo lunar seismic network was established during 
the lunar landing missions in the period 1969-72 and 
operated until it was shut down in 1977. Despite the very 
limited number of stations, the network was highly successful 
in determining the seismicity and internal structure of the 
Moon (Fig. 1) and in documenting the nature of inter- 
planetary objects in Earth-crossing orbits (e.g., Nakamura 
et ai, 1982). This success, however, did not come as a simple 
extension of what we knew about seismology on Earth. In 
fact, the majority of seismic events found on the Moon 
would not have been detected had we simply looked for 
normal earthquake-type signals with instruments of a sens- 
itivity typical for terrestrial use. 


SEISMIC STATIONS 



Fig. 1. Schematic diagram of the lunar interior inferred from the Apollo 
seismic data. The information depicted was derived from records obtained 
over an eight-year period by a four-station network. From Nakamura et 
al. (1982). 


The tidally triggered deep moonquakes, which constituted 
the majority of the detected events, occurred at depths 
significantly greater than those of any earthquakes. Rare 
but energetic shallow moonquakes, tectonic in origin, had 
a much higher frequency content than normal earthquakes 
of comparable magnitudes. All seismic signals, with long and 
nearly incoherent codas, were much more prolonged in signal 
duration than their counterparts on Earth (Fig. 2). More 
than 99% of the detected signals would have been too small 
to be detectable had we deployed instruments of sensitivity 
comparable to those used on Earth. That we were able 
to detect these events and use them quite profitably was 
because the instruments used on the Moon had sufficiently 
high sensitivity, at least two orders of magnitude higher 
than any on Earth, and because we collected continuous 
data without assuming a priori any characteristics of seismic 
signals or hypocentral locations for seismic sources. 

From this experience, it is clear that when we design 
and conduct a seismic experiment on a new planet such 
as Mars, we must (1) not assume that seismic signal char- 
acteristics are similar to those we see on the Earth or the 
Moon, (2) design an instrument with the highest possible 
sensitivity, and (3) be alert to the possibility of unexpected 
seismic sources. It is quite plausible that seismic signals carry- 
ing important information about the interior of Mars may 
come from such sources. In short, we must be prepared 
for the unexpected and be ready to capitalize on whatever 
signals are observed in such an experiment. 

Another important consideration in designing, a seismic 
network is the cost of operation. The Apollo network was 
terminated prematurely because it was too costly to keep 
operating, even though it was fully operational. With 
currently available technology, the operational cost of a 
martian seismic network can be reduced considerably from 
that of the Apollo network. Since a seismic network on 
Mars needs to be operational for an extended period of 
time because of the nature of the signals we aim to acquire, 
it is important that the routine network operation be made 
as simple as possible to minimize the cost of continued 
operation. 

THE VIKING SEISMIC EXPERIMENT 

The primary emphasis on Viking lander science was on 
biology, organic chemistry, imaging, and meteorology, with 
most aspects of surface chemistry, petrology, and geophysics 
relegated to future missions. The exceptions were the 
inorganic analysis experiment, the magnet experiment, and 
the seismometer. However, these studies were limited only 
to reconnaissance measurements. As a result, the Viking 
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Fig, 2. Representative Apollo lunar 
seismograms shown at a compressed 
timescale. From Nakamura et al. 
(1982). 
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seismic experiment was severely constrained by strict limi- 
tations on weight, power, and data rate, and was perturbed 
by the conflicting demands of the other onboard experi- 
ments. The weight constraint precluded an ultrasensitive 
seismometer of the Apollo class or a broadband seismometer. 
The original desire to place the seismometer on the surface 
was sacrificed because of the weight and complexity penalties 
of such an operation; thus an onboard location was dictated 
that immediately increased the noise level (by at least three 
orders of magnitude) because of lander and wind activity. 
The data-rate constraint required severe data compression 
using an onboard data processor and thus also imposed weight 
and power penalties. Overall, the most severe constraints 
on the Viking seismic experiment were the limited data 
allocations and the onboard location of the seismometer. 

These various constraints and trade-offs led to the de- 
sign of a short-period three-component seismometer with 
onboard data compaction and triggering to optimize the 
data return ( Anderson et al, 1972). The objectives were 
(1) to characterize the seismic noise environment at the 
landing sites, (2) to detect local events in the vicinity of 
the lander, and (3) to detect large events at teleseismic 
distances. Under optimal conditions it would also have been 
possible to determine the following: (1) the approximate 
distance of events from the separation of various seismic 
phases, (2) the direction of events to within a 180° ambiguity 
in azimuth, (3) the attenuation and scattering properties 
of the crust to determine if the crust were Moon-like or 
Earthlike in these characteristics (which are related to the 
volatile content), and (4) an estimate of crustal thickness 
if crustal and reflected phases could be identified. 


The Viking 1 seismometer failed to uncage, and no useful 
data were returned. The Viking 2 seismometer, emplaced 
on the surface of Mars in the Utopia Planitia region, 47.9°N, 
225.9°W, successfully uncaged and operated for a period 
of 19 months ( Anderson et al., 1977; Lazarewicz et al, 1981). 

The Viking experiment, even though it was highly 
curtailed by mission priorities and by the inability to uncage 
the first seismometer, provided valuable information for the 
planning of future missions to Mars; see also Appendix A. 
First, it established that the seismic background noise on 
Mars due to winds and atmospheric pressure fluctuations 
is very low. The Viking seismometer, mounted on a relatively 
compliant spacecraft with a large surface area exposed to 
winds, could still operate at maximum sensitivity at least 
half the time with no indication of noise visible on the 
records, a situation that indicated that seismometers with 
much greater sensitivities can be operated on the planet. 
Emplaced by penetrators or deployed as small packages, 
seismometers more sensitive than the Viking instrument 
by a factor of at least 10 4 can operate on the planet without 
being affected by typical martian winds. 

The indications from Viking are that Mars is probably 
less seismically active than the Earth (Anderson et al, 1977; 
Goins and Lazarewicz, 1979). Greater sensitivity is a must 
for any future seismic instruments on the planet. At the 
same time, a much broader frequency response is desirable 
than was used for the Viking seismometer, opening up the 
possibility of detecting normal modes, surface waves, tidal 
loading, and Chandler wobble. 

Another important consideration is the deployment of 
a network of instruments. With a well-placed network of 
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four very sensitive Apollo seismographs it was possible to 
study the seismicity of the Moon and to determine its 
internal structure, even though most moonquakes are very 
small. A similar approach with the deployment of a network 
of highly sensitive instruments is needed for Mars explo- 
ration. Since Mars is a larger planet than the Moon and 
shows much greater geologic and tectonic diversity, both 


global and regional seismic networks are needed to under- 
stand martian structure and tectonics. If Mars, as expected, 
is intermediate in seismic activity between the Moon and 
the Earth, then numerous seismic events of internal origin 
will be recorded by a martian seismic network. It is certainly 
not valid to say that the Viking experiment shows that 
Mars is an inactive planet. 


Expected Rates of Seismicity on Mars 
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The principal natural seismic events expected on Mars 
are marsquakes and meteoroid impacts. While the rates of 
occurrence are not well known for either type of event, 
theoretical considerations and experience from the Moon 
provide a basis for their estimation. We do not consider 
explicitly in this report possible man-made seismic sources 
on Mars, such as the impact of spacecraft onto the martian 
surface or the mechanical or explosive excitation of seismic 
waves by a lander or surface rover. The use of such sources 
would obviously augment an otherwise passive global seismic 
experiment. 

MARSQUAKES 

Volcanic and tectonic landforms on the surface of Mars 
provide evidence that the martian lithosphere has been a 
source of seismic activity over essentially the entire history 
of the planet. Images of young volcanic deposits that reveal 
only a few impact craters at the highest Viking Orbiter 
imaging resolution (e.g., Carr et al, 1977) suggest that 
internal activity may persist to the present. Thus, Mars 
should be a seismically active planet today, albeit at lower 
levels than in its more active geologic past. In this section 
we review several different potential sources of marsquake 
activity and present the results of calculations (see Appendix 
B) of the predicted rate level of seismic activity as a function 
of marsquake size. 

On a local or regional scale, potential sources of mars- 
quakes include magma motion (including extrusion), land- 
slides, time-varying insolation, and freeze-thaw cycles. Some 
of these sources will contribute to the seismic background 
noise. In other cases, characteristics of the seismicity may 
be diagnostic of the nature of the physical process involved 
(e.g., marsquake swarms near a volcano might be associated 
with magma motion). The largest of these events will be 
reliable seismic sources for studying interior structure. 

On a larger scale, marsquakes may result from the brittle 
failure of the lithosphere in response to stresses induced 
by such processes as crustal thickening by volcanism and 
plutonism, vertical motions arising from lithospheric buoy- 
ancy changes, and horizontal tractions exerted by convective 
flow in the mantle. Additional potential seismogenic pro- 
cesses on a large scale are the cooling and thickening of 
the lithosphere and changes in the principal moments of 
inertia. The rates of marsquake generation associated with 
this latter group of processes can be estimated (Appendix 
B) and provide a lower bound on the expected level of 
seismic activity. 

Differential cooling of the martian lithosphere constitutes 
an ongoing source of potential seismogenic strain. Such 
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Fig. 3. Expected rate of marsquakes generated by the global cooling of 
the martian lithosphere. The quantity N (M 0 ) is the number of marsquakes 
per Earth year with seismic moment equal to or greater than M„. The 
equivalent surface-wave magnitide Mg is also indicated, under the 
assumption the M s -M„ relation is similar to that on Earth. 

differential cooling is the major contributor to seismicity 
in young oceanic lithosphere on Earth ( Bratt et al, 1985; 
Bergman, 1986) and may be at least partly responsible for 
the generation of shallow moonquakes ( Nakamura et al., 
1979). For a simple parameterized convection model of the 
recent thermal history of Mars the rate of global cooling 
is equivalent to a mean lithospheric strain rate of 1.0 X 
10" 19 s' 1 (Appendix B). If the distribution of marsquakes 
by moment is similar to that in terrestrial oceanic lithosphere, 
then this mechanism yields a rate of seismicity as a function 
of marsquake size as given in Fig. 3. Global lithospheric 
cooling thus is expected to generate about 12 marsquakes 
per year with moment in excess of 1 0 25 dyn cm and 2.5 
events per year with moment in excess of 10 24 dyn cm. 
Events with moments in excess of about 10 23 dyn cm (about 
mt, 4.6) are routinely located on Earth with body-wave arrival 
time readings from global seismic networks, and for events 
with moments in excess of 5 X 10 23 to 10 24 dyn cm (about 
mj, 5) the source mechanisms are routinely determined (e.g., 
Dziewonski et al., 1987). 

Of course, other mechanisms will contribute to lithos- 
pheric strain and will thus add to the expected population 
of marsquakes over that depicted in Fig. 3. For instance, 
greater than average rates of cooling beneath the major 
volcanic provinces in Tharsis and Elysium will lead to 
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enhanced rates of marsquake activity in those areas. Solar 
tides and variations in polar cap loading due to changes 
in the planet’s obliquity contribute comparatively high strain 
rates over half cycles of l-ltP-yr duration, respectively, and 
may act to trigger marsquakes where the sense of strain 
reinforces that associated with longer-term processes (Ap- 
pendix B). Ongoing strains associated with mantle convec- 
tion or the motions of any fluids (magma or water) through 
the planetary interior will also act to raise the level of 
seismicity above that due solely to global cooling. 

We conclude that internal seismic activity will occur on 
Mars on a rate more than adequate to record, locate, and 
characterize marsquakes with a global seismic network. Over 
a 10-year period of observations, an effectively global 
distribution of activity is expected, although areas of 
significant topographic relief, large gravity anomalies, and 
comparatively recent volcanic activity are all favored sites 
for higher than average rates of seismicity. The records from 
such a population of marsquakes should be sufficient to 
address all the principal scientific objectives for a martian 
global seismic experiment described above. 


METEOROID IMPACTS 

Impacts of meteoroids in the mass range between 10 2 
and 10 6 g detected on the Moon by the Apollo lunar seismic 
network contributed significantly to our understanding of 
the deep internal structure of the Moon as well as the 
characteristics of the population of small objects in Earth- 
crossing orbits (Duennebier et ai, 1975). Similar observations 
on Mars are expected to be equally useful. 

There are several factors that influence the seismic 
detectability of meteoroid impacts. They include (1) the 
impact rate on the martian surface of various classes of 
meteoroids; (2) the effect of the atmosphere in retarding 
impacts; (3) the efficiency of impacts in generating seismic 
waves; and (4) the efficiency of seismic wave propagation 
through the interior of Mars. The combined effect of these 
factors in relation to the sensitivity of the sensors and the 
expected levels of ground noise, then, will determine the 
detectability of seismic signals from impacts. The effect of 
each of these factors is evaluated in Appendix C. 

Meteoroid impact rates depend on the population den- 
sity of various types of objects in the inner solar system, 
as well as their approach velocity and the size and mass 
of the target planet. Available data indicate that we can 
expect impacts of both cometary meteoroids and asteroidal 
fragments at rates, per unit area, on Mars similar to those 
on the Moon. Impact rates of asteroidal objects, however, 
may be somewhat higher on Mars than on the Moon if 
there is abundance of meteoroids originating from asteroids 
with Mars-crossing orbits. 


The effect of the martian atmosphere may be significant. 
Available data on this important effect are rather scarce, 
but they indicate that most meteoroids of cometary origin 
are likely to be effectively consumed in the atmosphere, 
while those of asteroidal origin in the mass range of interest 
are likely to be only partially affected. 

The efficiency of conversion of impact kinetic energy 
to seismic wave energy depends on the impact velocity, 
as well as other factors such as the properties of the target. 
The impact velocities of cometary objects are expected to 
be only slightly lower on Mars than on the Moon, while 
those of asteroidal origin are expected to be about 30% 
lower on Mars than on the Moon, or about half the kinetic 
energy for a given mass. The difference in physical properties 
between the martian surface and the lunar surface may 
be significant for seismic wave generation. However, we do 
not have reliable data at present to estimate this difference. 

There are substantial differences in seismic wave pro- 
pagation and attenuation between the Moon and the Earth. 
Available data indicate that the peak amplitude of the 
scattered seismic wave train observed on the Moon is about 
an order of magnitude greater than the amplitude of the 
P wave on Earth for a seismic source of a given energy, 
while the initial P-wave amplitude on the Moon is about 
an order of magnitude smaller than that on Earth. What 
it would be like on Mars is difficult to estimate, but the 
likely high volatile content of the martian interior suggests 
that seismic wave propagation in Mars may be similar to 
that in the Earth. As noted below, however, even if the 
intrinsic attenuation of seismic waves is similar in Mars and 
the Earth, the decay of seismic wave amplitudes with angular 
distance on Mars is less than on Earth because of the smaller 
planet size. 

As noted above, the Viking seismometer was not capable 
of determining the level of ground noise on Mars. However, 
it was also noted that when the wind was calm the observed 
ground noise was below the threshold of detection of the 
Viking seismometer. Thus, although we expect that wind 
stress will be the dominant source of noise, and that when 
the wind speed is high the ground noise will similarly be 
high, when the wind is calm the ground noise will be quite 
low. The ground noise at times of calm winds may be 
significantly below typical background noise on Earth, where 
seismic noise is dominated by wind-driven ocean waves. 

Given the greater surface area of Mars than the Moon, 
the likely higher seismic attenuation in the martian interior, 
the effect of the martian atmosphere, and the probable 
intervals of significant wind-generated ground noise, we 
expect a somewhat reduced detection level of large impacts 
and a greatly reduced detection of small cometary impacts 
on Mars compared with the Moon. An estimate of the 
rate of seismicity due to large impacts of asteroidal origin 
is derived in Appendix C and shown in Fig. 4. The rates 


Cumulative Annual impact Rate on Mars 
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Equivalent Body-Wave Magnitude 


indicated in the figure are conservative by as much as one 
to two orders of magnitude. Nonetheless, for a given size 
of seismic events, meteoroid impacts are expected to be 
less frequent than marsquakes. 


Fig. 4. Expected rate of meteoroid impacts per Earth year expressed 
in terms of equivalent body-wave magnitude mi,. The underlying 
assumptions are conservative, and the actual rates may be as much as 
a factor of 10-100 greater than those shown. 
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Seismic Observations in Relation to Objectives 


For each potential model for the internal structure of 
Mars, including assumed values for the radius and state of 
the core, the moment of inertia, the nature and location 
of any mantle discontinuities, and the thickness of the crust, 
seismic velocity distributions can be derived through the 
application of equations of state (Anderson, 1967, 1972) to 
a mineralogy satisfying the zero-pressure density. Okal and 
Anderson (1978) investigated both the properties common 
to all generalized models, as well as those that could resolve 
key unknowns, within the framework of seismological 
observations that could be made on the surface of the planet. 
Figures 5-7 present one example of such an inquiry, showing 
the computed density and seismic velocity profiles, the 
travel-time curves of the principal body waves traveling 
through the interior of the planet, and a sketch of the 
actual geometry of the direct P waves through the mantle 
and of PKP waves transiting through the core. In addition, 
Figs. 8 and 9 show the dispersion of Rayleigh and Love 
waves, which are dominantly at longer periods than body 
waves and travel along the surface of the planet. The periods 
of the various modes of free oscillation of the planet were 
also computed for a variety of models. 

MANTLE STRUCTURE 

Because of the smaller size and lower gravitational 
acceleration of Mars, the effect of pressure on density and 
elastic constants will be significantly weaker than on Earth, 
resulting in gentler gradients of seismic velocity. If the planet 
has a large temperature gradient (e.g., because a layered 
mantle prevents efficient convection), velocities could be 
practically constant throughout the mantle. If, on the other 



Fig. 5. Density and seismic velocity profiles for internal structure model 
AR. From Okal and Anderson (1978). 



Fig. 6. Travel-time curves for body wave phases P, PcP, PKP, S, ScS, 
and SKS in model AR. Note that PKP has been folded about the 180° 
axis. After Okal and Anderson (1978). 

hand, efficient convection has left the mantle at relatively 
cool temperatures, the mantle is olivine-rich, and the core 
smaller, phase transitions (similar to Earth’s 400-km and 
650-km discontinuities) could occur in the lower mantle, 
between 1000 and 2000 km depth. One such transition 
is shown in the model presented in Fig. 5, resulting in the 
triplication of P arrivals between distances of 50° and 90° 
(Figs. 6 and 7). 



Fig. 7. P and PKP ray paths inside Mars for model AR. A surface-focus 
event is assumed. The take-off angles vary from 1° to 50° in 1° increments 
and from 16.5° to 17.5° in 0.1° increments. The smaller diagram, 
reproduced from Julian and Anderson (1968; Fig. 9), is for the Earth. After 
Okal and Anderson (1978). 
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Mantle velocities may be obtained by inversion of observed 
travel times vs. distance by the Herglotz-Wiechert method, 
barring the presence of strongly developed low-velocity layers 
(see below). The resolution of mantle velocities, and in 
particular the recognition of possible stepwise increases, will 
be critical to an understanding of the composition, thermal 
state, and differentiation history of the planet.. It will be 
important in this respect to have good distance coverage, 
particularly in the distance range 50° -90°, to allow 
identification of possible triplications. 

SIZE AND STATE OF THE CORE 

While all models require that Mars have a core, the size 
and physical state of the core are poorly constrained: The 
core could be either liquid or solid, and the possibility of 
an Earth-like solid inner core cannot be discounted. If the 
(outer) core is liquid, S waves incident vertically on the 
core-mantle boundary (CMB) will be totally reflected, leading 
(as on Earth) to large-amplitude signals for core-reflected 
S (ScS) waves recorded at short distances (0°-20°) after 
bouncing on the CMB (Fig. 6). Also, in the case of a liquid 
core, the relatively low pressures involved lead to very low 
velocities in the core (4.6 km/s), a significant velocity drop 
at the CMB, and an extensive shadow zone for P phases 
between 100° and 140° (see Figs. 5 and 6), as well as for 
their S counterparts (direct S and SKS). 

If the core is solid, the higher P wave velocities in the 
core would substantially reduce the shadow zone for P and 
eliminate the one for S. In addition, the amplitude of ScS 
at short distances would be considerably Teduced. 

Crucial evidence for the state of the core will therefore 
be obtained from observations at close ranges (ScS, 0° -20° ) 
and at far distances (100° and beyond). In this respect, 
it is vital that the seismic network span a wide range of 
distances on the planet and, in particular, that some stations 
be located at considerable distances from areas of potentially 
greater than average rates of seismicity. 

If an inner core is present, significant chemical differ- 
entiation can be anticipated, with more of any lighter 
element (e.g., sulfur) in the liquid outer core. The inner 
core-outer core transition might then be sharper than in 
the Earth, and the resulting triplication would be detectable 
in the 90-180° distance range. 

Finally, the state of the core profoundly affects the 
fundamental modes of free oscillation of the planet, with 
the gravest spheroidal mode (qSj) having a period of about 
2100 ± 200 s if the core is liquid and only about 1350 
s if it is solid (Olcal and Anderson, 1978). Should sufficiently 
large marsquakes occur to excite these oscillations above 
ambient noise, their observation would strongly constrain 
the state of the core. 


CRUSTAL STRUCTURE AND VOLATILES 

The thickness of the martian crust can, in principle, be 
resolved by the inversion of first -arrival travel times at very 
short distances. Such an experiment would, however, require 
dense coverage along a profile of stations or sources, an 
unlikely geometry for a passive experiment given our inability 
to predict locations of activity. An alternative method is 
the use of surface wave (Rayleigh and Love) dispersion at 
periods in the range 10-100 s. The method can work both 
with the record at a single station of surface waves from 
a source at a known location and origin time ( Landisman 
et al, 1969) or over a small array of three or more stations, 
where the phase velocity can be measured directly across 
the array (e.g., T alandier and Okal, 1987) independently 
of the precise epicenter and origin time. In this respect, 
three-station broadband arrays would not only permit a 
direct measurement of the crustal structure, but their 
deployment at various sites could also help to identify lateral 
heterogeneity. 

Similarly, surface wave dispersion will be the key method 
to identify and resolve any low-velocity zone (LVZ) in the 
crust or upper mantle. LVZs cannot be formally obtained 
by the inversion of body wave travel times, so their resolution 
benefits significantly from the use of surface waves. The 
decrease in elastic moduli and increase in seismic attenuation 
evidenced by an LVZ (such as that between 100 and 250 km 
depth beneath oceanic areas on Earth) has usually been 
interpreted as due to partial melting, probably induced by 
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Fig. 8. Rayleigh wave phase and group velocities computed for model 
AR. 
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the presence of volatiles such as water. The mapping of 
LVZs, on a global or regional scale on Mars, would give 
insight into the history of planetary volatiles and could 
suggest links to magmatism and volcanic activity. 


Fig. 9. Love wave phase and group velocities computed for model AR. 
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The ability of a seismic station on Mars to detect a signal 
from a given seismic event depends on several factors 
independent of the intrinsic sensitivity of the sensors and 
recording system. These factors include the level of seismic 
noise, the rate of attenuation of seismic wave energy within 
Mars, and the decay of seismic amplitudes with recording 
distance. 

SOURCES OF SEISMIC NOISE 

Seismic noise levels on Mars are likely to be well below 
those detected on the Earth. Seismometers to be placed 
on Mars should therefore be designed to have the highest 
sensitivity possible to detect the smallest events possible. 
Some of the factors affecting relative levels of seismic noise 
on Mars, Earth, and the Moon, and instrument design con- 
siderations for reducing noise are given in Table 2. It is 
quite likely that even the most sensitive instruments and 
electronics (as are practical) will not be able to detect natural 
noise levels on Mars above their own intrinsic noise levels 
for much of the time and over a broad range of frequencies. 
If noise levels are as low as expected, it could be possible 
to detect very small marsquakes, possibly as low as mb 3, 
from anywhere on the planet. While natural noise sources 
such as wind stress and changes in diurnal thermal stresses 
will certainly contribute to the visible seismic noise at times, 
a major component of the expected noise could result from 
the presence of the equipment placed on the martian surface. 

Important natural noise sources are expected to include 
thermal cracking events caused by diurnal thermal stress 
variations in the near-surface rocks and soils, and noise 
caused by propagation of wind stresses into the martian 
surface. Thermal cracking events are seldom observed as 
seismic noise on Earth, mainly because other noise sources 
dominate and because thermal variations are relatively small 
(approximately 20 K/day, compared with about 100 K/day 
on the Moon and Mars). Thermal events were a major source 
of noise at frequencies above 3 Hz in the Apollo lunar 
seismic data. On Mars, such events will probably be even 
more common, as the length of the day is much shorter 
and stress variations will be more pronounced, although 
they will not penetrate as far into the surface. Major thermal- 
stress event swarms were observed on the Moon during 
eclipses, indicating that even very shallow penetration will 
cause these events. 

The lander and seismic package should be constructed 
to minimize thermal variations that could cause instrument- 
generated seismic noise. Insulation from direct sunlight 
should accomplish this. 


The high wind velocities on Mars will almost certainly 
contribute to the martian seismic noise levels. Wind 
impacting on the Viking lander was the only known external 
source of noise observed by the Viking seismometer. Wind- 
induced noise on Earth makes it necessary to avoid trees 
and other structures when emplacing seismometers on Earth, 
and burial in deep boreholes is preferred. While the wind 
velocities on Mars are often higher than those on Earth, 
the density of the atmosphere is much lower, and wind 
stresses will be less. However, motions of the lander and 
seismic package caused by the wind are still expected to 
be the largest single noise source. Thus, every effort should 
be made to minimize their cross sections and streamline 
their shapes. In addition, a substantial decrease in lander- 
generated noise can be realized by placing the sensor package 
at some distance from the station. 

TABLE 2. Seismic noise considerations on Mars. 

Instrument Design Considerations 


Detached from lander 
Cross section in atmosphere 
Burial 

Thermal insulation 

Earth Noise Lunar Noise 

(Compared to Mars) 

Atmosphere 



Density 

Higher 

NA 

Wind speeds 

Lower 

NA 

Surface roughness 

Same 

NA 

Q 

Same 

Higher 

Thermal Variations 

Lower 

Lower 

Volcanic/Seismic Activity 

Higher 

Lower 


SEISMIC ATTENUATION 

Seismologists have had the opportunity to measure seismic 
attenuation on two planetary bodies — the Earth and the 
Moon. The Earth is a tectonically active planet with an 
atmosphere and a hydrosphere, in contrast to the inactive, 
volatile-poor Moon. Attenuation is typically parameterized 
by the quality factor “Q.” High Q corresponds to low 
attenuation, and vice versa. The crust and upper mantle 
of the Moon have seismic Q values in the range 4000 to 
15,000 for shear waves in the frequency band 3 to 8 Hz 
(Nakamura and Koyama, 1982). For the Earth the corres- 
ponding values range from less than 25 to about 2000. At 
higher frequencies Q values in the Earth have been measured 
to 3000 and higher. Research on lunar materials indicates 
that the state of vacuum of the lunar environment is the 
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primary factor that dramatically increases Q on the Moon 
relative to the Earth (T ittmann et al, 1975, 1979). In this 
regard, the existence of an atmosphere and presumed sub- 
surface volatiles on Mars suggests that attenuation on Mars 
is more similar to the Earth than the Moon. 

The principal physical variables that influence Q for 
materials at depth (i.e., under pressure, with pores and cracks 
substantially closed) are temperature and pressure. At a 
constant pressure (depth), moderate temperature variations 
lead to strong Q variations since attenuation is an expo- 
nentially activated process. Thus, on Earth the hotter, 
tectonically younger provinces have systematically greater 
attenuation than older, more stable, and cooler regions. For 
example, a moderate-sized earthquake that might be felt 
in southern California would, for the same size event, be 
felt throughout the entire eastern United States if it had 
occurred there. Thus, higher temperatures and, indirectly, 
recency of tectonic activity lead to greater attenuation. 

To first order, then, the relative thermal environments 
of Mars and Earth likely dictate the relative attenuation. 
If radioactive elements are similarly distributed on Earth 
and Mars, then the thermal environment of Mars is likely 
to be similar to the stable shield provinces of Earth. The 
shield provinces on Earth are the low-attenuation, high- 
Q end members in the range of attenuation over the Earth. 
Under these assumptions, seismic wave propagation con- 
ditions on Mars are likely to be as good as the best conditions 
on Earth. If mantle heat transport has been sufficiently 
less efficient on Mars than on Earth so that mantle temper- 
atures are now higher, however, there may be greater atten- 
uation than beneath the Earth’s shields. 

For the core of Mars, the value of attenuation will depend 
strongly on its physical state. Liquid metals have nearly 
infinite Q values. Thus, a solid core on Mars implies a greater 
attenuation of seismic core phases as well as of free oscillation 
modes that are sensitive to core properties. 


AMPLITUDE-DISTANCE RELATIONS 

The rate of decay of amplitude with distance for seismic 
waves on Mars depends on several factors, including 
geometrical spreading and intrinsic attenuation; further 
details are given in Appendix D. Stronger body-wave 
amplitudes will be recorded on Mars than on Earth at a 
similar angular distance, A, from a seismic source of 
comparable moment, primarily because of the smaller size 
of Mars. Geometrical spreading alone accounts for a factor 
of 3.5; the effect of attenuation will depend significantly 
on the Q structure of the planet, which is still unknown. 
At a period of 1 s, if Q in Mars is at least as high as in 
the Earth, the final amplitude increase is a factor of at 
least 8.5 (and possibly as much as 15) for P waves and at 
least 50 (and as much as 3000) for S. The latter figure 
raises the possibility of a wealth of short-period teleseismic 
S data, in contrast to the case of the Earth, where large 
travel times and low Q values effectively eliminate them. 

For surface waves, geometrical spreading is in dependent 
of the radius of the planet and varies as 1 / \/sinA. Anelastic 
attenuation depends critically on the Q structure of the 
planet, but generally the shorter the distance traveled, the 
less the attenuation. On the other hand, surface-wave 
excitation is significantly larger on Mars than on Earth for 
a seismic source of a given moment and depth. As a whole, 
Rayleigh wave spectral amplitudes should be higher by a 
factor of about 3 (depending on Q) on Mars than Earth 
for a similar source. Time-domain amplitudes should be larger 
by an even greater factor, because of significantly less 
dispersion on Mars due to Jesser gradients of velocity with 
depth (see above). Love waves should be similarly enhanced 
on Mars for a given seismic source. 


Required Station Characteristics 
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SURFACE VS. BURIED SENSORS 

Technical specifications for the seismic instrument depend 
on the type of installation. Two installation modes are 
currently under consideration: surface installation by a rough 
lander (10-100 g decelerations) and subsurface installation 
using a penetrator. The pros and cons of these two modes 
are listed in Table 3. The trade-offs are as follows: 

Surface installations will be more susceptible to thermo- 
elastic strains. Given that the diurnal temperature variation 
is 100 K at the equator, the thermal inertia of the installation 
will buffer out short-term fluctuations, but long-term vari- 
ations (e.g., longer than 200 s) will appear as significant 
noise on a broadband instrument. Because of this effect, 
surface installations will probably be restricted to short period 
bands (0.04-20 Hz). 

A surface installation will be much more susceptible to 
wind-induced noise than a penetrator. This problem can 
be mitigated to some extent by removing the seismometer 
from the lander and placing it some distance away in a 
low-drag housing. A detached seismometer package has the 
added advantage of isolating the sensors from mechanical 
noise on the lander spacecraft, an option that would be 
difficult on a penetrator. Seismometer coupling will be much 
better with a penetrator than with a sensor emplaced on 
the surface unless some provision is made to attach a surface 
sensor with a spike or similar device. Burial of the seis- 
mometer is an ideal that may be met by penetrator instal- 
lation, but uncertainties in penetrability, equipment survival, 
and the technical difficulty in designing a power source 
(RTG/solar cell) that can withstand penetrator impacts (of 
order 10,000 g) are serious obstacles at this stage. Unless 
long-term power supplies tolerant to high g can be developed, 
or a mechanism of cushioning the power supply impact can 
be devised, penetrator installation with long-term power 
must be ruled out. 


TABLE 3. Comparison of surface and buried seismic instruments. 


Surface Installation 

Penetrator Installation 

pro: low-g impact (10-100 g) 
pro: relatively independent of site 
conditions 

con: thermoelastic strain 
con: short period only 
(25 s-20 Hz) 

pro: long-term power possible 

con: high g (10,000 g possible) 
con: site dependence critical- 
boulders, lava flows, alluvium 
pro: less sensitive to temperature 
pro: broad-band (DC-20 Hz) OK 

con: long-term power sources surviv- 
ing impact not possible with cur- 
rent technology 


TECHNICAL SPECIFICATIONS FOR A 
SURFACE SEISMIC STATION 

Specifications for an ideal surface seismic station on Mars 
are those obtainable in a first class observatory on Earth. 
These specifications are listed in Table 4. Were a penetrator 
installation feasible, specifications approaching an Earth 
observatory could be achievable for sites where the pene- 
tration was a few meters. If the schedule of station em- 
placement or limitations on station mass or complexity were 
overriding constraints, then — as discussed further below — 
some of the specifications could be relaxed for a subset 
of the network stations. 

On the basis of the discussion of sources of seismic noise 
given above, it is clear that surface sensors must be placed 
directly on the martian surface, i.e., off the lander. The 
greatest single improvement in signal-to-noise ratio is 
achieved by this action. Further improvements can be ob- 
tained by placing the sensors at some distance from the 
lander (preferably several spacecraft dimensions distant) and 
by burying the sensors within the soil (to lessen thermal 
perturbations). The need for distance between seismic 
sensors and spacecraft is heightened if the lander includes 
a meteorological boom, which will exacerbate the lander- 
generated seismic noise. 


TABLE 4. Seismic station specifications. 


Specification 

Ideal Seismic 
Station 

Martian 

Surface Installation 

Number of channels 

3 

3 

A/D sample rate 

1000 sps 

1000 sps 

A/D 

24 bit 

24 bit 

Stored data rate 

50 sps 

50 sps 

Data compression 

3:1 

3:1 

Data rate 

100 Mbits/day 

1 00 Mbits/day 

Bit error rate 

CIO' 8 

<10' 4 

Bandwidth 

DC - 30 Hz 

0.04 - 20 Hz 

Sensitivity 

Mass 

io-' 1 g 

io-‘°g 

Sensors only 

3 kg 

1.5 kg 

Sensors plus housing 

10 kg 

5 kg 

Volume 

0.02 m 3 

0.01 m 3 

Power 

2 W 

2 W 

Onboard RAM 

4 Mbytes 

4 Mbytes 

Calibration 

1/day 

l/day 
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DISCUSSION OF SPECIFICATIONS 

1. Number of channels. Three channels are required 
to discriminate longitudinal P-wave motion from transverse 
S-wave motion and to distinguish polarized surface waves. 
Teleseismic body waves have small angles of incidence so 
that teleseismic S is difficult to record well on vertical 
component records. 

2. A/D sample rate. A high sampling rate, 1000 sps, 
allows recursive anti-alias filtering, or other filters tailored 
to meet the site conditions, to be applied to the data stream. 

3. A/D. A 24-bit A/D provides 140 db of dynamic 
range, i.e., 10' 10 g to 10~ 3 g. 

4. Stored data rate. 50 sps/channel is typically used 
to record local earthquakes on Earth. 

5. Data transmission rate. This figure is calculated 
assuming continuous recording at 3 channels X 24 bits X 
50 sps X 1 day = 3.1 X 10 8 bits. It is estimated that with 
data compression this figure can be reduced by a factor 
of at least 3; see below. However, further compression, e.g,. 
by event detection, which in principle could reduce the 
total by a factor of 10 (to 10 7 bits), is much less attractive 
in that experience shows that array event detection on Earth 
often ends up with only half the array triggering. If all the 
data are collected, postacquisition processing can extract 
the signal from the remaining stations. 

6. Bit error rate (BER). Because of the finite 
bandwidth of the seismic sensor, a greater BER can be 
tolerated than for many other experiments. Therefore, a 
larger amount of high error-rate data is preferred to a small 
amount of error-free data. 

7. Banduridth. Die bandwidth depends on the ther- 
mal stability of the installation, how well coupled the sensor 
is to the surface, and how streamlined the superstructure 
is, in order to minimize wind effects. Critical to the per- 
formance of the site will be the lander-generated noise. If 
the lander includes booms that resonate in the wind at 
seismic frequencies, such motions will be the chief source 
of noise. Long-period waves may be difficult to observe at 
the surface. However, with an interrogate system and a 
broadband instrument, the bandpass might be altered if 
conditions (e.g., at night) are favorable, 

8. Power. Power is based on the requirement of the 
Guralp seismometers (low power version, 100 mW/channel) 
and the Reftek data acquisition package (the IRIS-designed 
portable digital seismograph) that operates at 1.5 W. Data 
acquisition systems with lower rates of power consumption 
(a few hundred mW) can be designed if power is a limitation. 

9. RAM. The RAM estimate is based on the assump- 
tion that data are transmitted to an orbiter every 8 hours, 


which for 10 8 bits/day gives a requirement of 4 Mbytes 
storage. 

10. Calibration. Calibration once per day would be 
used to track seismometer performance. 

DATA COMPRESSION 

The data compression scheme used by the Global 
Seismographic Network on Earth utilizes a first -difference 
encoding. With 24-bit data samples, most of the samples 
can be compressed to a one-byte first difference. When 
high-dynamic-range signals are encountered, either two or 
three bytes are used as necessary. The encoding uses bytes, 
rather than bits, for speed and convenience when using 
standard computers. The typical performance of this 
encoding is about 3 to 1 compression. These data may be 
further compressed by standard computer file-compression 
methods such as Huffman or Limpel-Ziv, achieving an 
additional 33% reduction. Thus, seismic data compression 
of about a factor of 5 is possible. 

In considering compression schemes, the robustness of 
the methods in the presence of bit errors must be considered. 
The standard first -difference encoding method used by the 
terrestrial Global Seismographic Network packages data into 
64-byte frames, each of which contains an absolute reference 
sample. Thus, if there is a bit error, no more than about 
60 bytes of data are lost. With alternative coding schemes 
that do not “packetize” the information, bit errors can 
introduce substantial data loss in compressed data, and this 
consideration is important in the final selection of a 
compression scheme. 

STATION COMMAND ABILITY 

Because we do not know the form that all seismic signals 
or noise will take on Mars, it will be extremely important 
to have the capability to send a few simple commands to 
each seismic station after acquisition and analysis of initial 
data. Commands can include control of corner frequencies 
and rolloff rates of bandpass filters, the application of notch 
filters or other noise suppression algorithms, changes in 
sampling rates or instrument gain settings, or modifications 
to data compression parameters. As the flexibility of the 
instrument is increased and the adaptability to unforeseen 
seismic signals is maximized, trade offs are possible in the 
required overall rate of data transfer from each station to 
an orbiter relay. With sufficient station commandability it 
is possible that an overall data rate of 10 Mbits/day per 
station might be sufficient to achieve all the scientific 
objectives of a global seismic network on Mars. 


Station Siting Requirements 
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THE STATION TRIAD CONCEPT 

A station siting plan must strike a compromise among 
the various scientific objectives. There must be a mix of 
closely spaced stations to detect and locate nearby seismic 
events and to serve as arrays for measurements of phase 
velocity and propagation direction, and more globally 
dispersed stations to provide planetary coverage and to 
ensure the recording of core phases and other signals 
diagnostic of deep structure. While broadband, three- 
component systems are preferred for all stations, we can 
envision that, because of different emplacement systems 
(e.g., penetrator vs. lander), cost constraints, or sequential 
station emplacement scenarios, there may be strong reasons 
to consider a simpler and more rugged station with only 
short-period seismometers. The recommended plan, con- 
sistent with overall mission constraints as currently defined, 
consists of nested triads of stations. 

Small triads should consist of stations located approx- 
imately 100 km apart in local arrays. If there must be a 
compromise among station types, then we recommend that 
each small triad consist of 3 three-component short-period 
instruments arranged in a triangle around a broadband 
observatory-class three-component sensor. If broadband 
three-component sensors can be emplaced at all sites, then 
the central station in each small triad could be omitted. 
If this second scenario is followed, in order to provide backup 
signals in the event of an unexpectedly hard landing or 
high levels of long-period noise, we recommend that the 
broadband sensors be augmented at each station with an 
additional set of more rugged short-period seismometers. 

Three of these small triads should be emplaced in a larger 
triangular pattern, approximately 3500 km on a side. One 
such large triad should occupy each of the eastern and 
western hemispheres of Mars. This arrangement gives a total 
of 12 stations in each hemisphere, or 24 stations in all. 

RATIONALE FOR STATION SPACING 

The 100-km spacing of instruments in the small triads 
will permit the detection (by comparison of signal character 
on several stations) and location (by inversion of local P 
and S travel times) of seismic events within and near the 
triangular area, so several of the small triads should be located 
in regions likely to have above-average rates of seismic 
activity. For instance, the Viking 2 seismometer was 
emplaced about 110° from the Tharsis region, possibly the 
most seismically active region on Mars. It has been estimated 
that an instrument 20 times more sensitive could have 
detected at least 120 seismic events per year if it were located 


closer to Tharsis (Goins and Lazarewicz, 1979). The station 
spacing for the small triads represents a compromise between 
resolution of source location and areal coverage. A simple 
rule of thumb in local networks is that the depth of focus 
of a natural seismic source can be resolved as long as the 
nearest station is closer to the epicenter than one focal 
depth. Most earthquakes are located in the upper 20 km, 
but the lesser thermal gradients expected on Mars than 
on Earth and the likelihood that thermal stress associated 
with global cooling will generate marsquakes throughout 
an elastic lithosphere as much as several hundred kilometers 
thick suggest that a 100-km spacing will be adequate to 
locate most of the local events recorded by each small triad. 

In addition, a 100-km spacing will permit the triads to 
serve as quadripartite arrays for measurement of surface- 
wave phase velocities, which for periods of 15-100 s have 
wavelengths of 50-450 km (Figs. 8 and 9). The triads will 
also serve to measure the propagation azimuth and phase 
velocity of body wave phases. Among the key body wave 
phases at short distances that will benefit from phase velocity 
measurements are near-vertical reflections (which arrive at 
very high phase velocities) from the crust-mantle boundary, 
from the core, and potentially from any prominent mantle 
discontinuities, as well as direct and mantle-refracted waves 
from local events (which tend to arrive at phase velocities 
closer to the seismic wave speed) that yield information 
on crustal structure. Azimuth and phase velocity measure- 
ments of waves from events distant from all triads will 
constitute primary information for the location of such 
sources. 

The 3500-km spacing of triads in each hemisphere will 
permit the siting of a number of triads in distinct regions 
likely to have greater than average rates of seismicity. The 
3500-km distance, about 60° of arc, will ensure that some 
stations are located at a key distance range to record body 
wave triplications from mantle discontinuities. In addition, 
the three triads will ensure a reasonable level of detectability 
and locatability of seismic sources throughout each hemi- 
sphere. 

The siting of one large triangle in each hemisphere will 
give effectively global coverage. Further, such an arrange- 
ment will afford the opportunity to observe core phases 
from seismic sources well located on the basis of travel times 
to one or more of the triads on the opposite hemisphere. 

STRAWMAN STATION LOCATIONS 

Our suggested strawman plan for seismic station layout 
includes two hemispheric arrays — one in the eastern hemi- 
sphere and one in the western hemisphere (Fig. 10). All 
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Fig. 10. Possible locations (triangles) for seismic station triads for a martian global seismic network. The background map depicts topography (in 
kilometers) with respect to a fourth degree and order equipotential. 


locations are at altitudes less than 6 km, as required by 
present mission constraints, and all but one are at absolute 
latitudes less than 45° to permit the use of solar power. 
The proposed western array has vertices in Tharsis west 
of Ascraeus Mons (15°N, 1 10° W), just north of the eastern 
portion of Valles Marineris (7.5°S, 52.5°W), and in the 
southern highlands (45°S, 105°W). The choice of Tharsis 
is obvious, since it is the region of the planet with arguably 
the most recent volcanic activity ( Tanaka et al., 1988) and 
consequently highest interior cooling rates. Valles Marineris 
has high topographic gradients and displays evidence for 
recent volcanism ( Lucchitta , 1987) and landslides. A small 
triad in the southern highlands extends western hemisphere 
coverage to the south (e.g., activity triggered by time-variable 
polar loads) and provides a key locale for the determination 
of regional structure (e.g., crustal thickness) to test models 
for the origin of the crustal dichotomy. 

The proposed eastern array has a vertex in the Elysium 
province (22°N, 211°W), in the Isidis Basin (5°N, 270° W), 
and in the northern lowlands (65°N, 275° W). Elysium is 
a major volcanic province displaying evidence for relatively 
recent activity ( Plescia , 1990). Isidis is a mascon basin (Sjogren, 
1979); shallow moonquakes have shown some tendency to 
be associated with mascon mare basins on the Moon 
(Nakamura et al., 1979). The placement of a triad in the 
northern lowlands should permit the measurement of local 


crustal structure, making the site complementary to that 
in the southern highlands for tests of the crustal dichotomy, 
and will extend coverage to the northern portions of the 
planet, including the north polar region. 

STATION SEQUENCING SCENARIOS 

Mission design considerations may dictate that stations 
in a martian seismic network be emplaced in a sequence 
spanning several distinct launches and perhaps several years. 
The most important constraint for system emplacement is 
that as much of the array as possible must record data 
simultaneously. This constraint dictates that design lifetimes 
of individual stations must be considerably greater than the 
time interval over which the network is established. 

If seismic activity on Mars consists of numerous large 
events, and the first stations emplaced are widely scattered 
(e.g., several thousand kilometers apart), then observations 
from a sparse, early network could be used to locate regions 
on the planet where seismic activity is greatest and where 
small triads emplaced later in the sequence could be located 
for optimum results. If, in contrast, seismic activity on Mars 
occurs mainly as small, well-distributed events poorly 
recorded at distant stations, then a sparse network of widely 
spaced stations may be able to locate few, if any, events. 
This second situation is the more likely and favors the early 


EAST 
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emplacement of a small triad of stations in a single area 
(e.g., Tharsis). Observations from such an early array could 
yield the location of near events, determine local structure, 
and assess the size distribution and approximate locations 
of distant large events. Because we do not know the character 
of martian seismic signals, it is essential that two or more 
of the earliest emplaced seismic stations be located near 
(i.e., within 100 km of) each other so that signals from 
tectonic and impact events can be recognized and distin- 
guished from wind-generated and thermal noise. 

We note one additional advantage to a sequenced station 
emplacement scenario in which signals from the first stations 
are received and analyzed in time to affect the operational 


characteristics of later-emplaced stations. Data rate is likely 
to be a strong constraint on experiment design, particularly 
once the full network is operational. While many data 
compression schemes can be envisioned for seismic data, 
as discussed above, implementation of the more efficient 
of such schemes is extremely hazardous in the absence of 
information about the characteristic durations, bandwidths, 
and frequency of martian seismic events. We expect that 
the analysis of signals from the first seismic stations on Mars 
will permit implementation of data compression algorithms 
that will considerably reduce the necessary data rate from 
the final, full network with little or no loss of scientific 
return. 
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RECOMMENDATIONS FROM THE 
VIKING SEISMOLOGY TEAM 

(from Lazarewicz et ai, 1981) 

The Viking Seismology Experiment 

The purpose of the Viking Seismology Experiment was 
to determine the seismicity of Mars and define its internal 
structure by detecting vibrations generated by marsquakes 
and meteoroid impacts. After obtaining more than 2100 
hours (89 days) worth of data during quiet periods at rates 
of one sample per second or higher, the Viking 2 seismo- 
meter was turned off as a consequence of a lander system 
failure. During the periods when adequate data were 
obtained, one event of possible seismic or meteoroid impact 
origin was recognized; however, there is a significant 
probability that this event was generated by a wind gust. 
The lessons learned from Viking, however, will ensure that 
seismic systems on future Mars missions will have a con- 
siderably higher probability of obtaining their goals. 

The Viking seismometer was a three-component short- 
period system designed to meet severe constraints of weight, 
size, power consumption, and data rates necessary for 
incorporation into the Viking lander (Anderson et ai, 1972). 
Its small size and location on the lander body produced 
a relatively insensitive and noisy seismic system. The Viking 
seismometer is about 1/20 the sensitivity of the Apollo 
seismometers at 3 Hz. Over most of the frequency band 
of seismological interest, the Apollo seismometers are 
generally 2 to 3 orders of magnitude more sensitive than 
the Viking seismometer. Because of lander vibrations in 
response to martian winds, many of the seismic data are 
contaminated and thus unsuitable for detection of seismic 
events. When the winds are quiet (<2 m/s), background 
noise is below the level detectable by the system. An event 
of possible seismic origin (SOL 80 event) was recognized 
at a time when wind data suggest relatively quiet conditions. 
(One SOL is one martian solar day = 24 hr 39 min 35.25 s.) 
No wind measurements were made within 20 minutes of 
the event, so it is possible that a gust' could have occurred 
between wind samples and produced the event ( Anderson 
et ai, 1977). If the event is seismic in origin, the high 
frequency and short duration of this event suggest that 
it was generated locally. 

The detection of one local marsquake (at most) during 
the equivalent of 89 days of operations allows us to set 
limits on the probable seismicity of the planet. The prob- 
ability is greater than 67% that Mars has a lower seismicity 
than the Earth’s intraplate seismicity. If martian seismicity 


is similar to intraplate seismicity on Earth (there is no 
evidence of plate tectonics on Mars), then about 2 to 3 
events would have been detected (Goins and Lazarewicz, 
1979). 

Requirements for successful seismic experiments on Mars 
and other planetary objects can be met using present-day 
technology. These include (1) high sensitivity, dynamic 
range, and frequency bandwidth to allow detection of small 
and distant events; (2) seismometer networks to allow 
location of detected events and inversion of seismic param- 
eters to obtain planetary structure; and (3) flexible data 
collection and compression methods to allow variation of 
parameters for optimum retrieval of information. 

Much like the Viking seismometer on Mars, the first 
seismometer on the Moon (Apollo 11) was noisy and told 
us little about the Moon. A great wealth of information 
was obtained, however, by later Apollo seismometers, and 
we are confident that the same will be true when more 
advanced seismic systems are installed on Mars. 

In evaluating the Viking Seismology Experiment, it is 
important to know the history of the experiment, and thus 
understand the reasons for the design limitations. 

Viking was designed primarily as a biologically oriented 
mission. The search for extraterrestrial life was the moti- 
vating force, and other scientific goals were secondary. The 
original tentative instrument list did not include a seismo- 
meter. While design changes were taking place for the Viking 
lander, a seismometer was proposed and accepted for the 
Viking mission. From the outset, the primary constraints 
for this instrument were weight, power, and cost, in that 
order. Initial tests, using the Caltech engineering unit and 
a crude lander model, indicated no significant difference 
among placing the seismometer on, under, or near the lander 
for seismic coupling of the instrument to the surface. Any 
seismometer deployment mechanism would have exceeded 
the weight allocation for the seismic experiment and would 
have forced a decision between a lander-mounted seismome- 
ter or no seismometer at all. Through testing and com- 
promise, the original design evolved into the instrument 
that went to Mars as a 2.2-kg, 3.5-W instrument. The weight, 
power, data allocation, and cost constraints forced a 
relatively low-sensitivity, narrow-band, survey experiment. 
Although the “noise” environment of Mars is known to 
some extent, we do not know the characteristics of internal 
events or the effect of propagation on these characteristics. 
The high-Q surface scattering layer in the Moon results 
in very emergent (signals having a gradual onset) seismic 
arrivals, long codas that obscure secondary arrivals, and the 
destruction of coherently dispersed wave trains (Goins, 1978). 
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The lunar seismic signatures have proven to be substantially 
(and unpredictably) different from the Earth’s, leading to 
the obvious inference that seismic signatures need not be 
similar from one planetary object to another. 

The Apollo seismometer systems collected signals that 
were so radically different from the Earth’s that they were 
not identified as seismic signals until after the second Apollo 
(12) mission; even then identification was possible only due 
to the high data rates (40 X 10 6 bits/Earth day) and effective 
real-time interaction between scientists and instruments. 
A similar case may be made for Mars. The fact that there 
were no positively identified signatures does not imply that 
there are no seismic events on Mars. The Viking seismo- 
meter, if placed on the Moon instead of the Apollo seismic 
systems, would have detected lunar seismic activity at a 
rate of only 0.5 event per year, compared to approximately 
20,000 events per year observed with Apollo instrumen- 
tation. In other words, the Viking instrument probably 
wouldn’t have detected any seismic activity on the Moon. 

Recommendations for Future Seismic 
Experiments on Mars 

The Viking Seismology Experiment was a valuable source 
of experience in performing extraterrestrial seismology. This 
was the first seismic experiment on an unmanned, extra- 
terrestrial mission. It was expected that the constraints on 
mass, power, and data allocation and the physical location 
on the spacecraft would compromise the scientific results. 
The amount of compromise and specific problems became 
well understood only after the landings on Mars. In this 
sense, one of the major accomplishments of this experiment 
was an improved ability to design similar future missions. 
Four specific recommendations were presented by the Viking 
Seismology Team. All four are technically feasible and will 
significantly improve future extraterrestrial seismic experi- 
mentation. 

1. Isolation from noise. We know of three noise 
sources on Mars: lander activity, thermal noise, and wind 
noise. A seismic package should be isolated from all three 
as well as possible. 

This is the most important recommendation. A first 
seismic experiment on a planetary object searches for seismic 
signals of unknown character. Any detected signal must 
be thoroughly analyzed to determine its origin and thus 
classify it as a seismic signal or noise. A landed spacecraft 
is necessarily noisy, and the analysis and categorization of 
the noise adds greatly to the effort, time, and money spent 
for data analysis. 

The emplacement of a seismometer on a planetary object 
is not as straightforward a problem as on Earth. Construction 
of a seismic vault is presumably impossible; even searching 
for a good site (if one can get to it) is a difficult process. 


It is reasonable to assume the existence of a surface layer, 
probably of unknown or poorly understood properties, sepa- 
rating the surface from bedrock (if any). This surface layer 
will probably modify, possibly severely, an incident seismic 
signal (as it did on the Moon), thus making interpretation 
of seismic signals that much more difficult. Furthermore, 
the surface layer undergoes diurnal changes of temperature 
with the associated thermal fluxes and changes in the 
thermal state of a surface seismometer. 

Constructing a seismometer package showing a small cross- 
sectional area to the wind, and planting it on the surface 
of Mars, should increase the threshold for wind noise to 
300 m/s. The sensitivity of the seismometer could be 
increased by 4 orders of magnitude and the seismometer 
would still be free of wind noise ( Anderson et al, 1977). 

The ultimate goal is total isolation from spacecraft and 
surface noise sources. In order of importance, a planetary 
seismometer must have its sensors (a) physically separated 
from the spacecraft, (b) buried, and (c) attached to bedrock 
(if any). It may be possible to couple seismometer emplace- 
ment with a coring mission. After a core is removed from 
a hole dug by a lander, the hole could be filled with an 
instrument package containing seismic and other sensors. 
An instrument package could also be buried with a corer. 
Alternatively, a penetrator mission could easily carry a 
seismometer and bury the instrument at each landing site. 

2. Seismometer networks in seismically active 
areas. The Earth and the Moon have been found to have 
seismic and aseismic zones. Mars is also expected to be 
heterogeneous in seismic activity. Most notably, the Tharsis 
area is expected to be characterized by stresses that should 
produce seismic activity even if the area is in isostatic balance 
( Sleep and Phillips, 1979), so seismic activity in this area 
is likely to be higher than elsewhere on the planet. There 
are no planetary objects where a uniform surface seismicity 
is expected. The strength of seismic signals received by a 
seismometer generally decreases with distance; as a result, 
a strong bias exists for sensing the seismicity of the immediate 
area. For reasons of spacecraft landing safety, the Viking 2 
lander was placed in Utopia Planitia, where substantial 
seismic activity was not expected. Theoretical considerations 
of martian seismicity, and the physical limitations of the 
Viking seismometer, place 80% of the potentially detectable 
events within 10° (590 km) of the lander. It has been 
estimated that an event in the Tharsis region would have 
to have had a magnitude greater than 9 to be detected 
by the Viking 2 seismic experiment (Goins and Lazarewicz, 
1979). Obviously, the location of a seismometer on a plan- 
etary object is very important. 

A single seismometer, especially without detailed know- 
ledge of the nature of seismic signals, is insufficient for deter- 
mining the location and nature of the seismic source. A 
network of seismometers offers many advantages over a single 
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instrument, depending on how the network is deployed. 
Three useful types of seismic networks are 

(a) Subwavelength network: Sensors closely spaced 
(within a fraction of a seismic wavelength) so that 
the time of arrival of a signal at the different instru- 
ments is nearly simultaneous. The signals may be 
added, reinforcing the coherent signals through con- 
structive interference, while weakening the incoherent 
signals (noise). In this fashion, the signal-to-noise ratio 
is improved as the square root of n, where n is the 
number of sensors. 

(b) Local network: Sensors are widely spaced, 
relative to a seismic wavelength, so that the arrival 
time differences are much greater than the measurable 
time resolution, but the instrument separation is much 
less than the planetary radius. With proper adjustment 
of the relative phases of incoming signals, this phased 
network becomes a “steerable beam” for distant 
events, where azimuthal direction can be calculated 
for a given signal, or a preferential azimuth for a par- 
ticular study can be chosen. Local seismic events can 
be located directly by triangulation and/or travel-time 
difference analysis. 

(c) Planetary network: The network is spread 
throughout the entire planetary surface. This network 
allows for studies of heterogeneity in seismicity and 
interior structure, as well as focal mechanism studies. 
This approach has the added feature (and compli- 
cation) of having instrumentation in different geo- 
logical provinces. 

By far, the best system is a combination of all three types. 
If practical constraints do not allow for a combined system, 
some type of network is nevertheless desirable. The simplest 
network is one seismometer with every landed package, a 
c-type array (such as Apollo). A penetrator mission could 
deploy a planet-wide c-type network. A planetary rover could 
emplace a b-type network, for example, three seismometers 
separated by 1-10 km. Obviously, the more members in 
the network, the better the resulting resolution. 

3. Flexibility through software control. The ability 
to control instrumental parameters from Earth greatly in- 
creases the flexibility of the whole system and could be 
used to filter in (or out) particularly desirable or undesirable 
signals. 

The sampling rate of an instrument depends on the 
frequency band of interest and data constraints. The samp- 
ling rate must be more than twice the upper frequency 
limit to prevent aliasing. The ability to modify the recorded 
sampling rate could be used to save data space. If data 
compression is used, similar software modification of the 
degree of compression could also be used. 

4. Event detection, data collection, and compres- 
sion. An important aspect of seismic data is that more 


than 90% of the data is relatively uninteresting. Important 
information comes in bursts at random and unpredictable 
times. While Apollo could afford to send all its seismic data 
back to Earth, where relative importance could be deter- 
mined at the leisure of scientists, Viking did not (and future 
planetary seismic experiments probably will not) have this 
luxury. Decisions on information content will have to be 
made before the data are returned to Earth. It is important 
to maximize the scientific content in the available data 
space; however, it will probably not become evident how 
best to accomplish this goal until after the experiment 
package is placed on another planetary object. 

In order to sample three axes at a sampling rate of 20 Hz, 
8 bits per sample over a 24-hour period, a seismometer system 
requires 41.5 X 10 6 bits plus overhead (e.g., data manage- 
ment, time codes). Estimating 3% of the data space for 
overhead, or 1.2 X 10 6 bits, a total of 42.7 X 10 6 bits are 
required. The Viking seismometer was allocated a nominal 
1.1 X 10 6 bits per SOL, in other words, a 39:1 reduction 
from optimum in data allocation. There are three imme- 
diately obvious approaches to this problem: (1) increase the 
data allocation, (2) devise a data compression technique 
where only the most important parts of a given seismic 
event are kept, and (3) devise a priority system of data 
storage that eliminates unwanted noise and periods of low 
information when more informative data are encountered. 

The total data allocation for a given spacecraft is entirely 
constrained by the state of the art of data storage and 
transmission capabilities. This is a spacecraft system con- 
straint and is mostly a problem of engineering and economics. 
Dividing the total data allocation among the various experi- 
ments will, for the foreseeable future, leave any seismology 
experiment far short of the optimum 42.7 X 10 6 bits/day. 
This forces optimization of the in situ data collection, 
processing, compression, and transmission schemes. 

This recommendation, although a very important one, 
has no known straightforward or obvious solution. Many 
different approaches are possible, and much research must 
be done well ahead of the next mission so that a system 
may be developed and ready when the time for hardware 
development arrives (normally around the time of mission 
approval). It is important to note that the primary target 
for this study is in situ data management. 

Conclusion 

The most severe problems encountered in the Viking 
seismic experiment were, in order of importance: (1) failure 
of the Lander 1 seismometer; (2) wind noise due to the 
location of the seismometer on the spacecraft; (3) inadequate 
data compression techniques; and (4) low sensitivity. 

It is a straightforward task to eliminate all four problems 
for future extraterrestrial seismometry experiments. The 
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Viking Seismology Team has put forward two sets of recom- 
mendations for future work on Mars: 

(1) Scientific goals - to determine 

(a) seismicity 

(b) velocity-depth relationships 

(c) nature of the core 

(d) location of partially molten areas 

(e) tectonic state of the Tharsis area 


(2) Engineering goals - to achieve 

(a) isolation from noise 

(b) emplacement of seismometer networks 

(c) flexibility in software control 

(d) improvement of in situ data handling 

All these recommendations are necessary and technically 
feasible. 
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Appendix B 


EXPECTED RATES OF MARSQUAKES 

Roger J. Phillips 


Introduction 

In this appendix we review several regional and global- 
scale processes that are potential sources of martian seis- 
micity, and we present the results of a few simple calculations 
to predict the level of seismic activity as a function of mars- 
quake moment and magnitude. 


Moment Recurrence Relation 

The most fundamental measure of marsquake size, and 
that most readily relatable to physical processes, is the scalar 
seismic moment M 0 given by 

M 0 = /iAu (1) 

where /a is the shear modulus in the marsquake source region, 
A is the fault area that ruptures during a marsquake, and 
u is the average coseismic fault slip. The distribution of 
seismic events by moment is generally given by a relation 
of the form 


N(M 0 ) = aM 0 - b (2) 

where N is the number of events of moment greater than 
M 0 , and a and b are empirical constants. We assume that 
equation (2) holds as well for marsquakes, and we further 
assume that the value of b (but not that of a) is identical 
to that for earthquakes in oceanic intraplate regions, a 
terrestrial tectonic setting possibly analogous to that of the 
martian lithosphere. From a catalog of oceanic intraplate 
earthquakes of known moment over the inclusive years 
1977-1988 (E. A. Bergman, personal communication, 1991), 
the b value is 0.67, and we adopt this figure for Mars. 

From equation (2), the number of events falling in the 
interval M 0 to AM 0 + M 0 is given by 

N (M 0 , M 0 + AM 0 ) = a [M 0 ’ b - (M 0 + AM o r b ] (3) 

The sum of the seismic moments of all events occurring 
over a time interval t can be expressed as 

M max 

[ 0 dN = f ° M 0 [a b M^ 1 *] dM 0 

J° ( 4) 

= ^ b W- b 



where M™ ax is the largest seismic moment for Mars. For 
a given mechanism of lithospheric strain, the characteristic 
rate of shear strain e can be related to the sum of the 
seismic moments over the time interval t through the 
relation ( Bratt el aL, 1985) 

eV — 2M 0 /pit (5) 

where V is the seismogenic volume. Thus 

± = i=±±(M, rf-'eptV 

t b (6) 

and the frequency of occurrence of marsquakes in a given 
interval of moment is given by the relationship 

N (M 0 , M 0 + AM 0 ) = a[M 0 - b - (M 0 + AM o r b ] (7) 

With an assumed value of /a, an estimation of strain rate 
and seismogenic volume yields a prediction of the frequency- 
moment distribution of marsquakes. This can be converted 
into a magnitude recurrence rule using a moment-magni- 
tude relation (Appendix D). 


Thermoelastic Mechanisms 

We first consider global cooling of a planetary lithosphere 
as a source of strain rate. If is planetary radius and L s 
is the thickness of the seismogenic lithosphere, then 


• R at f Ro 2 , 

t = * — I r dr 

Ro Rts/zLs -V 


( 8 ) 


where Ris the rate of change of R<„ a is the volumentric 
coefficient of thermal expansion, Tis the characteristic rate 
of cooling of the lithosphere, and R^ is the radius at mid- 
depth in the seismic lithosphere. 

A parameterized convection calculation has been used 
to predict the change in the martian lithospheric temper- 
ature profile during the last 100 Ma. This model includes 
the effect of the decay of radiogenic heat sources in the 
crust as well as a decrease with time in heat flux from 
the sublithospheric mantle. On the basis of analogy with 
the terrestrial oceanic lithosphere, we define the base of 
the seismogenic lithosphere at the 800°C (1073 K) isotherm 
(Wiens and Stein, 1983; Bergman, 1986). In the parameterized 
convection solution, this isotherm corresponds to a depth 
at present of approximately 150 km. The average temper- 
ature change within the seismogenic lithosphere during the 
last 100 Ma is about 1 1 K, and this figure is used to estimate T 
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The principal change in the temperature gradient with time 
is caused by the thickening of the lithosphere. For values 
of a and n of 3 X 10' 5 K" 1 and 7 X 10 10 Pa, respectively, 
the characteristic lithospheric strain rate is -1.0 X 10’ 19 s' 1 . 
Table B1 shows the predicted number of marsquakes per 
year as a function of seismic moment interval; a maximum 
moment of 10 27 dyn cm is assumed on the basis of the 
data for terrestrial oceanic intraplate earthquakes. The 
cumulative number of marsquakes per year is shown in Fig. 3. 


TABLE Bl. Predicted number of marsquakes per Earth year as a function 
of moment interval (dyn cm) due to global cooling. 


Logio Moment Interval 

Number/yr (Predicted) 

21-22 

200 

22-23 

42.8 

23-24 

9.15 

24-25 

1.96 

25-26 

0.42 

26-27 

0.09 


Thus, from the thermoelastic cooling of the martian 
lithosphere alone, a seismic network of 10-year duration 
could potentially measure about 116 events between 
moments of 10 23 and 10 27 dyn cm. Table B2 expresses this 
information in terms of surface wave magnitude Mg under 
the assumption that the M s -M 0 relation is that given for 
the Earth in Appendix D 

M s = log 10 M o - 19.46 (9) 

where M 0 is in dyn cm. From this relation, the maximum 
moment corresponds to an M s of about 7. As noted in 
Appendix D, the surface wave magnitude as measured from 
time-domain amplitudes may actually be somewhat greater 
for a given seismic moment on Mars than on Earth, so 
Taible B2 is conservative. For a 10-year network, there are 
predicted to be about 60 events between Mg 4 and 7 due 
to global cooling. 


TABLE B2. Predicted number of marsquakes per year as function 
of surface wave magnitude M$ due to global cooling. 


Magnitude Interval 

Number/yr (Predicted) 

1-2 

461 

2-3 

98.5 

3-4 

21.0 

4-5 

4.50 

5-6 

0.96 

6-7 

0.21 


The rate of lithospheric cooling, and thus the volumetric 
rate of marsquake generation, may be greater in the Tharsis 
area than for the planet as a whole. We have made a simple 
estimate of this effect on the supposition that the magmatic 


source region beneath Tharsis has been cooling for at least 
several hundred million years. A simple one-dimensional 
diffusion model is used. The magma source region is presumed 
to lie beneath a lithosphere of thickness L, which tracks 
the cooling beneath with a linear temperature gradient. As 
the source region begins to cool, the lithosphere is anom- 
alously thin; it thickens with time toward its ambient (global 
average) value. The change in temperature gradient is con- 
trolled by the rate of thickening of the lithosphere, which 
is given by 



where Tl is the temperature at the base of the lithosphere, 
T 0 is the anomalous temperature (Tl + AT) in the sub- 
lithosphere at t = 0, T s is the surface temperature, and k 
is thermal diffusivity. The background or ambient tem- 
perature gradient is given by dT amb /dz and is time dependent. 
Equation (10) has been used to evaluate the strain rate 
at mid-depth in the seismogenic lithosphere beneath Tharsis; 
the results are shown in Table B3. 


TABLE B3. Strain rate and thickness of seismogenic lithosphere 
due to cooling of a sublithospheric temperature anomaly, 

AT, of 300 K beneath Tharsis. 


Time (Ma) 

Strain rate (s ') 

Seismic Lithosphere 
Thickness (km) 

100 

-7.8 X 10' 19 

99 

200 

-3.6 X 10~ 19 

111 

300 

-2.4 X 10‘ 19 

116 

400 

-1.9 X 10~ 19 

120 

500 

-1.6 X 10" 19 

123 

600 

-1.5 XIO' 19 

125 

700 

-1.4 X 10‘ 19 

127 

800 

-1.3 X 10' 19 

129 

900 

-1.2 X 10' 19 

130 

1000 

-1.2 X 10“ 19 

131 

Ambient 

-1.0 x 10' 19 

154 


Thus, for several hundred million years after initial cooling 
of a 300-K thermal anomaly, the strain rate and thus the 
volumetric rate of seismicity in the vicinity of Tharsis could 
be more than twice that due to whole-planet cooling. 

Loss of Buoyancy Support 

Another process specifically associated with Tharsis 
involves loss of buoyancy support of the lithosphere. If the 
region beneath the elastic lithosphere is supporting the 
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Tharsis load in part by thermal buoyancy, then as this region 
cools, a load develops, which must adjust isostatically or 
flexurally. The density change in the source region is given 
by 


where G is the gravitational constant, M s is the mass of 
the Sun, and r s is the Mars-Sun distance, given by 

r s = a (1 - e cos E) (17) 


Ap = -<*PoT (11) 

where p 0 is the ambient density. The corresponding strain 
rate is given by 


• 6 A P g K _ a Po g h s T (12) 

EE E 

where h s is the thickness of the cooling source region and 
E is Young’s modulus. For a source region 400 km thick 
that has cooled 200 K in 500 Ma (Table B3), the resulting 
strain rate is about 10 20 s' 1 , a value small compared with 
the effect of lithospheric thickening. 


where a and e are the semimajor axis and eccentricity of 
the orbit, respectively. The term E is the eccentric anomaly 
and can be approximated by 


E = Eo = M + esinMH sin 2M 

2 


with 


(18) 


M = ^ L (t-T) 


(19) 


where P is the period of the orbit and T is the time of 
perihelion passage of Mars. The radial tidal displacement 
is given by 


Change in Principal Moments 

The cooling and contraction of the lithosphere in the 
Tharsis region will also change the difference in the principal 
moments of inertia, C-A, and this change will lead to 
membrane strain in the elastic lithosphere. 

For a spherical cap of half-angle y, density p, and thickness 
h, which is displaced downward by x 

A(C- A) = 47t x P h cos y (cos 2 y-1) (13) 


<5r = 


hW 

go 


( 20 ) 


where go is the martian gravitational acceleration, and h 
is the tidal Love number. The strain rate in Mars associated 
with this process is given by 


-3 G h 

go r s 


M s Ro 


8r s 

dt 


( 21 ) 


where Rq is the radius of Mars. The change in planetary 
flattening is given by 

3 A(C-A) 

Af = 

2 M (14) 

where M is the mass of Mars. On the equator 

e =“0.4 Af = 0.6— — Ro x P h cos y (cos 2 y-1) (15) 

M 

For h = 200 km and AT = 10 K, x 20 m, and with y = 45° 
and a 100 Ma time interval, e«-5x 10‘ 23 s' 1 , so this effect 
appears to be negligible. 

Annual Solar Tide 

Due to the eccentricity of the martian orbit, the solar 
tide is time variable with an annual period. The solar tide- 
raising potential is 

GMsR ° ( 16 ) 


and dr s /dt is obtained from equations (17-19). The resulting 
strain rates, of order 10' 17 s’ 1 , are high compared with the 
other processes that have been thus far examined. However, 
as discussed below, the total strains over one half cycle 
are small, and solar tidal effects are seen to be generally 
unimportant. 

Obliquity Changes and Polar Cap Loading 

The final mechanism examined for global lithospheric 
strain is that of polar cap loading due to the 10 5 -year obliquity 
cycle (Ward, 1979). Rubincam (1990) describes "post-glacial 
rebound” on Mars due to the changing obliquity and the 
consequent growth and decay of large COj loads on the 
martian poles. As a result of the polar cap changes, the 
difference in the principal moments (C - A) periodically 
changes as well, although there is a phase lag due to the 
finite viscosity of the mantle. This gives rise to a variation 
in the second harmonic of gravity, J 2 , with a perturbation 
magnitude given by 

M™ 

{11) 
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where is the mass of the C0 2 load and 'P is the half- 
angle of the cap deposits. The strain resulting from this 
process is given by 

e~0.4 Af = 0.4y AJ 2 (23) 

and the strain rate is given by 



where T is the period of the obliquity variation. R ubincam 
(1990) estimates a total load of M a p = 10* 7 kg with a half- 
angle 'P of 10°. With the 120,000-year obliquity variation, 
this gives a strain rate of about 1.5 X 10“ 19 s' 1 . 

Conclusions 

Several geophysical processes leading to a wide distribu- 
tion of strain rates have been examined. Those processes 
associated with thickening of the martian lithosphere can 
yield significant levels of seismic activity over a decade-long 
observing period. Other processes, such as periodic loading 
by obliquity variations, lead to comparable or higher strain 
rates. However, the total strain should also be considered 
before concluding which processes might contribute most 
to seismic activity. If we consider the strength of the 
lithosphere to lie between 10 and 100 MPa, then with a 
Young’s modulus of E of 10* 1 Pa, a total strain of 10" 3 
to 10 4 is required for marsquakes. The strain rates for the 
various mechanisms considered can each be multiplied by 
a characteristic duration of the process to estimate the total 
strain that might accumulate. These results are given in 
Table B4, and it is seen that only those processes associated 


with lithospheric cooling can, in and of themselves, lead 
to sufficient strain levels for marsquake generation. Several 
of the processes with higher strain rates (e.g., solar tides) 
may nonetheless act to trigger marsquakes if strain levels 
due to other processes are favorable. Overall, the mech- 
anisms considered here suggest that a martian global seismic 
experiment is quite feasible and that the expected level 
of marsquake activity is sufficient to achieve the principal 
scientific objectives of the experiment. 


TABLE B4. Summary of strain rates, durations, and total strains for 
the various mechanisms considered for marsquake generation. 


Process 

Strain Rate (s ') 

Duration of Pro- 
cess (Earth yr) 

Total Strain 

Thermoelastic 
stress: Whole- 
planet cooling 

1 X 10‘ 19 

1 X 10 9 

3.2 X 10‘ 3 

Tharsis: Cooling 
(AT = 300 K, 
t = 400 Ma) 

2X 10' 19 

4X 10 8 

2.5 X 10‘ 3 

Tharsis: Loss of 
buoyancy support 
due to cooling 

1 X ur 20 

4 X 10 s 

1.3 X 10 -4 

Tharsis: Change 
in principal 
moments 

5 X 10' 23 

1 x 10 9 

1.6 X 10-* 

Annual solar tide 
effect 

1 X 10 17 

1 x io° 

3.2 X 10-‘° 

Obliquity varia- 
tion: Polar cap 
loading 

2 X 10~ 19 

5X 10 4 

3.2 x 10 -7 


Appendix C 
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SEISMIC DETECTION OF METEOROID 
IMPACTS ON MARS 

Yosio Nakamura 

Meteoroids, those numerous small bodies in interplanetary 
space too small to be detectable by astronomical observa- 
tions, collide with planets and satellites that happen to be 
on their orbital paths. They range in mass from less than 
10 6 g, called micrometeoroids or interplanetary dust, to 
more than 10 15 g, thus approaching the size of asteroids. 
Those in the middle of this mass range, generally 10 2 to 
10 6 g, are of particular interest to seismologists. Smaller ones, 
though detectable by a seismometer if they fall close to 
the instrument, do not generate sufficiently large seismic 
waves to be useful in studies of regional to global scale, 
while larger ones are too rare to be expected within a 
reasonable time of observation (D uennebier et al, 1975). 

Impact signals of meteoroids are useful in two different 
ways. First, as these seismic signals propagate though the 
deep interior of a planet, they provide important information 
on the structure of a planetary interior (e.g., Nakamura, 
1983). Second, the spatial and temporal distribution of 
impacts on a planetary surface provides information on the 
orbital distributions of these interplanetary objects (e.g., 
Oberst and Nakamura, 1987). 

For these impacts to be useful in seismic studies on Mars, 
we must detect a significant number of impact seismic signals 
during the operation of the seismic stations. Can we 
reasonably expect to do so? To answer this question, we 
need to examine several factors, all of which influence the 
seismic detectability of meteoroid impacts. They include ( 1 ) 
the impact rate, or the rate of encounter with Mars, of 
various classes of meteoroids; (2) the effect of the atmosphere 
in retarding impacts; (3) the seismic efficiency, or the 
efficiency of conversion of impact kinetic energy to seismic 
waves; and (4) the efficiency of seismic wave propagation 
through the martian interior. The combined effect of these 
factors, then, must be compared with the sensitivity of the 
seismic instruments and the expected level of ground noise. 

There is a large quantity of information relevant to an 
assessment of these factors. They include ( 1 ) terrestrial obser- 
vations of meteors and meteorite falls; (2) astronomical 
observations of comets and asteroids; (3) the Apollo lunar 
seismic observations, which provided data on impact rates, 
temporal and spatial distribution of impacts, seismic 
efficiency, and seismic wave propagation for the Moon; 
(4) lunar and planetary crater statistics; (5) lunar chrono- 
logy, which assigns absolute ages to various lunar terrains 


of known crater density; and (6) theoretical results on the 
orbital dynamics of interplanetary objects. 

In the following we examine, from available information, 
these factors influencing the seismic detection of meteoroid 
impacts on Mars. Particular emphasis is placed on a 
comparison with the Moon, for which extensive information 
on the seismic effect of impacts is available. 

Impact Rate 

To estimate the rate of meteoroid impacts on Mars, we 
need to know what families of objects presently exist in 
orbits crossing that of Mars. A generally accepted, but not 
completely understood, scenario is that most meteoroids in 
the inner solar system are derived from collisional break- 
up and perturbation through encounters with planets of 
two principal types of objects, namely comets and asteroids. 
What proportion of these objects is cometary is difficult 
to determine from currently available observational data 
because all observational methods are strongly selective in 
detecting different kinds of objects. Because of this difficulty, 
impact rate estimates based on one type of observations, 
e.g., terrestrial fireball observations, may not be directly 
applicable to the estimation of impact rates relevant to 
seismic effects. 

Probably the most closely relevant set of data for our 
purposes comes from the Apollo lunar seismic network. 
During the 8 years of network operation, 1743 seismic events 
clearly identified as meteoroid impacts were detected by 
the long-period instruments ( Nakamura et al., 1982). [The 
short-period instruments detected several times more impact 
events (Duenrtebier and Sutton, 1974), but these events have 
not been fully analyzed.] A recent study by Oberst and 
Nakamura (1991) shows that many of the relatively less 
energetic impacts, too weak to be detected by multiple 
stations and believed to be of cometary origin, fall in streams, 
while more energetic impacts, believed to be either asteroidal 
or possibly derived from short-period comets, are mostly 
sporadic, with a few important exceptions occurring in 
streams or swarms, indicating relatively recent breakup from 
their parent bodies. Of these 1743 impacts, only 18 were 
large enough to be useful for investigation of the internal 
structure of the Moon by conventional methods (Nakamura, 
1983). 

To estimate the impact rate on Mars by inference from 
the lunar data, we need to know how these objects are 
distributed in the inner solar system. A Monte Carlo cal- 
culation of the fate of objects from a wide class of starting 
orbits by Wetherill (1975) shows that they evolve to pro- 
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duce a similar impact flux on all the terrestrial planets. 
This is in agreement with the observation that the cratering 
records on the Moon, Mercury, and Mars are similar. The 
actual impact flux, however, depends on the orbits of the 
source objects. Those in original Mars-crossing orbits, for 
example, produce about 10 times more impacts on Mars 
than on the Moon ( Wetherill , 1975). Thus, it is difficult 
to make a precise estimate of the ratio of impact fluxes 
on Mars and the Moon without knowing the sources of 
these meteoroids. On the other hand, observations of actual 
impact rates on both the Moon and Mars will give us clues 
as to the origin of these objects. 

Crater statistics on the Moon and terrestrial planets (e.g., 
Hartmann, 1973, 1977; Neukum and Wise, 1976; Neukum 
and Hiller, 1981: BVSP, 1981) provide a vast quantity of 
data that may be relevant to this problem. Unfortunately, 
only for the Moon do we have an absolute chronology of 
impact terrain through radioisotope dating of samples. Thus, 
the crater statistics by themselves do not provide an 
independent means of estimating the impact rate on Mars. 

For a given population density, n p , and approach velocity, 
Voo, of a class of meteoroids, the impact rate, F p , depends 
on the size of the planet on which the impacts occur 

F p = n p Voo S 2 /4R 2 

where R is the geometric radius of the planet, S is the 
effective radius given by 

S 2 = R 2 [l + (8 ttCRV3voc 2 )] 

and p is the mean density of the planet. 

The approach velocity v„ depends on the orbits of the 
planet and the impacting object. Table Cl lists estimated 
values of some relevant quantities for the Earth, Moon, 
and Mars as adopted from those given by Hartmann (1977) 
for typical asteroidal and cometary objects. Also listed are 
the corresponding impact velocities, v„ given by 

Vi 2 = vj + v e J 

where v esc is the escape velocity. 


TABLE Cl. Estimated values of approach velocity, effective planetary 
cross section, and impact velocity for the Earth, Moon, and Mars. 




Asteroids 



Comets* 


Voo 

(km/s) 

S 2 /4R 2 

<5 

(km/s) 

koo 

(km/s) 

S 2 /4R 2 

V\ 

(km/s) 

Earth 

14 

0.42 

18 

38 

0.27 

40 

Moon 

14 

0.26 

14 

38 

0.25 

38 

Mars 

8.6 

0.33 

10 

31 

0.26 

31 


Includes short -period comets; assumes v„ = 0.4 J'orb- 
* Long-period comers; assumes Va> = 1.3 


Atmospheric Effects 

The martian atmosphere, even though tenuous compared 
with that of the Earth, is sufficient to cause significant 
ablation and deceleration to meteoroids entering it. Gault 
and Baldwin (1970) have estimated that for a typical entry 
speed of 10 to 15 km/s, appropriate for asteroidal objects, 
those with masses less than about 10 3 g are decelerated 
and consumed in the upper atmosphere, and only those 
with masses greater than about 10 7 g are relatively 
unaffected. They estimate that the latter objects would 
produce craters larger than approximately 50 m in diameter. 
As mentioned above, most of the objects detected by the 
Apollo lunar seismic network are in the mass range from 
10 2 to 10 6 g. Thus, most of them would be expected to 
be strongly affected by the atmosphere of Mars. 

Since ablation increases significantly with increasing speed, 
objects coming from highly eccentric cometary orbits will 
be more effectively consumed by the atmosphere. Cometary 
objects are also likely to be more friable and thus more 
easily disintegrated in the atmosphere. The relatively low 
density of cometary material also contributes to higher 
ablation through the atmosphere because of the larger sur- 
face area for a given mass. Although exact calculations are 
not available, we expect that among those in long-period 
cometary orbits only rare, massive fragments and some 
higher-density objects sometimes found in cometary showers 
(Halliday, 1988) have any chance of being detected by a 
martian seismic network. 


Seismic Efficiency 

Only a very small fraction of the kinetic energy carried 
by the impacting object is converted to seismic energy, with 
most of the energy spent in excavating an impact crater 
and heating the target material. For example, the impact 
of spent spacecraft on the lunar surface during the Apollo 
missions produced an estimated total radiated seismic energy 
about 10 -6 of their preimpact kinetic energy ( Latham et 
al, 1970). The efficiency of the energy conversion depends 
upon many factors such as mass, speed, and density of the 
impacting object, angle of impact, and physical properties 
of the target material (e.g., Schultz and Gault, 1975). 

The diameter, D, of impact craters generally follows a 
power-law scaling with respect to the kinetic energy, Ej, 
of the object 

D = CE> 

where the exponent A is about 1/3. [An experimental value 
by Gault (1973) gives A =0.370]. If we assume that the 
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amplitudes of seismic waves scale similarly, then the 20% 
to 30% reduction in the expected impact speed from the 
Moon to Mars (Table Cl) translates to about a 15% to 
20% reduction in seismic amplitude, probably insignificant 
relative to other uncertainties. 

The density of the impacting object does not appear to 
affect the seismic efficiency greatly for the Moon. This is 
evidenced by the fact that both low-density cometary and 
high-density asteroidal objects are well detected by the lunar 
seismic network ( Oberst and Nakamura, 1989). An empirical 
relationship by Gault (1973) also shows only a very weak 
dependence of crater diameter on the density of the impactor 
with an exponent in a poweT-law scaling of 1/6. As stated 
above, the martian atmosphere is expected to have a much 
more profound effect on low-density objects than this slight 
variation in seismic efficiency. 

The physical property of the target material may play 
a significant role in determining the seismic efficiency of 
impacts. However, we do not presently have reliable data 
to estimate the difference in physical properties between 
the Moon and Mars. 

Seismic Wave Propagation and Attenuation 

Whether an impact of a given kinetic energy is seismically 
detectable at a given distance from an impact site on a 
planetary surface also depends strongly on how seismic waves 
propagate through the interior of the planet. There is a 
significant difference in the mode of seismic wave prop- 
agation between the Earth and the Moon. On the Earth, 
generally distinct P and S body-wave arrivals are followed 
by often dispersive surface-wave trains. On the Moon, in 
contrast, a relatively weak initial body- wave arrival is 
followed by a long train of scattered waves of much larger 
amplitude and little coherency. This unusual behavior of 
seismic wave propagation in the Moon is attributed to the 
existence of a highly heterogeneous near-surface zone of 
very low seismic attenuation. 

The observed amplitude decay with distance for P waves 
on the Earth and that of scattered body waves on the Moon, 
however, are remarkably similar. This probably reflects the 
fact that both the Earth and the Moon have a distinct 
low-velocity crust overlying higher-velocity mantle, thus 
producing similar geometrical spreading of seismic waves on 
both bodies. The amplitude generally decays rapidly with 
distance for the first 20° to 30° of distance, and then decays 
only slightly from 30° to 90°, beyond which more rapid 
decay follows (Veith and Clawson, 1972; Nakamura et al, 
1976). 

An empirical surface-focus amplitude-distance curve 
based on large explosions (Veith and Clawson, 1972) gives 
an average value for log(A/T) on Earth of about 1.5 for 


distances between 30° and 90°, where A is the P-wave 
amplitude in nanometers normalized to a source of body- 
wave magnitude m b = 5, and T is the wave period in seconds. 
m b = 5 corresponds to about a 100-kT yield (Bache, 1982), 
or 4.2 X 10 14 J of released energy. In contrast, the impacts 
of spent S-IVB rockets, each of which had a preimpact 
kinetic energy of about 4.6 X 10 10 J, produced peak ground 
vibration on the Moon of about 10 nm at 0.5 Hz at distances 
beyond 30°. Thus, normalizing for the same energy level 
as above, we obtain log(A/T) = 2.7 for peak amplitudes 
expected for impacts on the Moon. This means that the 
peak amplitudes of scattered wave trains on the Moon are 
about an order of magnitude greater than those of P-wave 
arrivals from surface explosions on Earth. The initial P- 
wave arrivals on the Moon, however, are about two orders 
of magnitude smaller than the maximum amplitudes of 
scattered wave trains. Therefore, they are about an order 
of magnitude smaller than those on Earth. 

The amplitude-distance relations for Mars are difficult 
to estimate without knowing the scattering and attenuation 
characteristics of the martian interior. If Mars is like the 
Earth, i.e., relatively low in scattering and low in Q, then 
amplitude decay should be similar to that on Earth. If Mars 
is like the Moon, i.e., intensive scattering and extremely 
high Q, then we may expect an amplitude decay similar 
to that on the Moon. However, if we have intensive 
scattering, which might be expected because of the past 
impact history, and lower Q, which might be expected 
because of a significant volatile content of the interior, then 
the decay of amplitudes with distance may be more drastic 
on Mars than on either the Earth or the Moon. 

Ground Noise 

For lunar seismic observations we were in a sense very 
fortunate that we had a combination of extremely high- 
sensitivity seismometers and extremely low ground noise, 
except for short intervals just after sunrise and sunset. (The 
Apollo seismographs had a peak sensitivity of better than 
0.1 nm in ground displacement at 0.45 Hz.) As a conse- 
quence, we detected literally hundreds of impact events 
every year. Whether we will see a similar result on Mars 
is not certain. 

The results of the Viking seismic experiment showed that 
when the wind speed was high, the seismic noise level was 
also high. However, the Viking seismometer was mounted 
on the lander, not on the ground, and its sensitivity was 
about two orders of magnitude lower than that of the Apollo 
lunar instruments. Therefore, there are no experimental 
data on how low the ground noise may be when the local 
atmosphere is calm. It may be possible, however, to estimate 
this quantity from available data. 
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Expected Rate of Seismicity from Impacts 

From the above discussion, it is apparent that most of 
the impact signals large enough to be useful for seismic 
investigation of the martian interior will come from spor- 
adic impacts of relatively large meteoroids of asteroidal origin. 
The frequency of occurrence of such impacts can be 
estimated from lunar data. 

The latest estimate of the flux of large sporadic impacts 
on the Moon (J. Oberst and Y. Nakamura, personal com- 
munication, 1990) gives an annual rate N of impacts of 
preimpact kinetic energy greater than E (in J) over the entire 
lunar surface of 

logioN = -0.99 log 10 E+ 11.4 

for 2 X 10 u < E < 2 X 10 12 J (mass range from 2 X 10 6 to 
2 X 10 7 g using an impact velocity on the Moon of 14 km/s, 
Table Cl). If we assume that the population density of 
large meteoroids in Mars-crossing orbits is the same as those 
in Earth-crossing orbits, the impact rate on the martian 
surface, after adjustment for the differences in impact 
velocities and effective radii as given in Table Cl and in 
the surface areas of the Moon and Mars, will be given by 

logjoN = -0.99 logioE + 11.6 

If we assume an efficiency of seismic energy conversion 
of 10” 6 as observed for the impacts of spent spacecraft on 
the lunar surface ( Latham et al., 1970) and the relationship 


logioE = 2.3 mt - 0.5 between energy and body-wave mag- 
nitude mb ( Stacey , 1977), we obtain an estimate of the rate 
of impacts in terms of equivalent magnitude as 

logioN = -2.3 mb + 6.1 

This relation is shown in Fig. 4. This estimate is likely 
to be conservative. Seismic efficiency observed for missile 
impacts at White Sands, New Mexico (R. Eggleton, as quoted 
by Latham et al., 1970) was an order of magnitude larger 
than that observed for the Moon. Furthermore, the pop- 
ulation density of large meteoroids in Mars-crossing orbits 
may be significantly higher than assumed here. Thus, the 
rate of seismogenic impacts may be significantly higher (by 
a factor of as much as 10 to 100) than shown in this diagram. 

Conclusions 

Although definitive data are still lacking, given the greater 
surface area of Mars than the Moon, but also the likely 
higher seismic attenuation in the interior, the effect of the 
martian atmosphere, and the probable intervals of significant 
wind-generated ground noise, we expect a somewhat reduced 
detection of large impacts, mainly of asteroidal origin, and 
a greatly reduced detection of small, cometary impacts on 
Mars compared with the Moon. 


Appendix D 
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EXCITATION OF SEISMIC WAVES 
ON MARS 

E mile A Ok al 

Despite the fact that the internal structure of Mars is 
not known, several general conclusions regarding the exci- 
tation of body and surface waves by seismic sources can 
be drawn from what is known and from terrestrial ex- 
perience. As a basis for simple calculations we use model 
AR of O kal and Anderson (1978) for the internal structure 
of Mars, and we assume that seismic sources consist primarily 
of double couples (i.e., marsquakes). 

BODY WAVES 

We refer to Kanamori and Stewart (1976) for the follow- 
ing expression of the amplitude of a body wave recorded 
at teleseismic distance A from a source 

u = M 0 — ' — G(A)CRA(t,Q) (1) 

4zrpv J 

We will define and discuss each of these factors separately. 

M 0 is the seismic moment, scaling the seismic source. 
Since the purpose of this computation is to assess the 
response of the planet to a unit moment source, we assume 
that the range of seismic moments is common to both 
planets. 

The next factor, ^7, where p is the density and v the 
seismic velocity (either P or S) at the source, characterizes 
the efficiency of the source material at generating body 
waves. This term should not be significantly different in 
Earth and Mars. 

The next factor represents the geometrical spreading of 
the body wave in a radially inhomogeneous, spherical planet 
as a function of angular distance A. For surficial sources 
G can be written (Kanamori and Stewart, 1976) 


The prominent difference between the two planets comes 
from the smaller size of Mars, through the a~ 3/2 factor in 
equation (2). The relative behavior of the other terms in 
equation (2) is more difficult to assess. Mars is probably 
a less heterogeneous planet than Earth (because of the 
reduced compression of the mantle), and crustal velocities 
can be assumed to be comparable. Thus, at the same angular 
distance, A, the rays are less strongly bent, i.e., sample a 
shallower fraction of the planet’s radius than in the Earth, 
resulting in stronger values of the incidence angle i, and 
therefore of sin i/cos 2 i. On the other hand, because of the 
lesser refraction and of the shorter travel times, the second 
derivative d 2 T/dA 2 is less than in Earth (by about a factor 
of 1.5). Given this complexity, we computed numerically 
G (A) in the case of P waves, both for model AR in Mars 
and for the Jeffreys-Bullen model in the Earth, and found 
that at typical mantle distances (25° to 95°), the martian 
value of G is about 3.5 times larger than the terrestrial 
one. Because neither planet is expected to feature drastic 
variations in Poisson’s ratio throughout its mantle, the ray 
paths for P and S waves will be nearly identical and the 
same result would be expected for S waves. 

The next two factors in equation (1) are C, which is 
the surface response to the incident wave, expected to be 
equivalent in Mars and the Earth, and R, the radiation 
coefficient characterizing the geometry of the double couple 
with respect to the surface and the azimuth of the station, 
also expected to be equivalent to existing geometries on 
Earth. 

The last factor in equation (1) characterizes the anelastic 
attenuation along the path of the body wave 


A (t, Q) — exp 
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where a is the planetary radius, and i is the take-off (and 
emergence) angle of the ray at the surface: i = sin' 1 pv/a 
with p = dT /dA, where T is the travel time. 


where to is the angular frequency, v(r) and Q(r) are the 
velocity and quality factor as functions of radius, and the 
integration is conducted along the wave path. Obviously, 
the values of A in Mars will depend critically on the Q struc- 
ture of the planet, especially for high-frequency waves. Even 
if Mars had the same Q as Earth, the shorter travel times 
T in equation (3) would boost A significantly, since the 
average mantle P wave on Mars travels about 420 s, as 
opposed to 600 s on Earth. Using an average Qp of 600 
for P waves at 1-s period, A would jump from about 0.043 
on Earth to 0.11 on Mars, a factor of 2.5. In the case 
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of S waves (Qs Earth = 250), 1-s waves would be attenuated 
50 times less than on Mars, simply on account of the shorter 
travel times. In addition, we might expect a somewhat higher 
average Q value on Mars than on Earth because of a less- 
developed asthenosphere. Using an average Mars Qp of 800 
and Qs of 450, we find A = 0.192 for P waves (a 4.5-fold 
increase with respect to the Earth), and 3000 times less 
attenuation for 1-s S waves. Finally, for “long-period” body 
waves (T = 10 s), the difference in attenuation is negligible 
(1.16) for P, and about a factor of 2 for S. 

In conclusion, and due primarily to various effects re- 
sulting from the smaller size of the planet, we find that, 
at a similar angular distance A from a seismic source of 
comparable moment, stronger body wave amplitudes will 
be recorded on Mars than on Earth. Geometrical spreading 
alone accounts for a factor of 3.5; the effect of atten- 
uation will depend significantly on the Q structure of the 
planet. At a period of 1 s, the final amplitude increase 
is a factor of at least 8.5 (and possibly as much as 15) 
for P waves; at least 50 (and as much as 3000) for S. The 
latter figure raises the possibility of a wealth of short-period 
teleseismic S data, in contrast to the case of the Earth, 
where large travel times and low Q values effectively elim- 
inate them. 

SURFACE WAVES 

Using model AR, we have computed excitation coef- 
ficients for spheroidal and torsional normal modes in the 
period range 15-430 s, and for source depth h ranging 
from 10 to 100 1cm. 


Rayleigh Waves 

We refer to Fig. 8 for a plot of the phase and group 
velocities of Rayleigh waves. An interesting characteristic 
of this diagram is the fact that the maximum dispersion 
occurs around a period T = 50 s, while inverse dispersion 
takes place for 20-s waves. This reflects the particular choice 
of 50-km crustal thickness in model AR. As a result, the 
most prominent component of a Rayleigh wave in the time 
domain may be centered around 50 s, rather than 20 s, 
as documented in the synthetic seismogram on Fig. Dl. 
Also evident is the absence of the inversely dispersed branch 
of mantle modes, with the longer periods actually arriving 
early in the record. 

Following Okal and Talandier (1989), we study the loga- 
rithmic excitation L R of Rayleigh waves averaged over many 
geometries of source and receiver. Fig. D2 shows the 
logarithmic average excitation for 10 depths ranging from 
10 to 100 km. This figure shows that the value for h = 
30 km can be used to model all selected depths with an 
error of less than 0.2 units of magnitude above T = 45 
s. A cubic spline fit to the curve for this depth yields a 
formula of the type 

L R = - 1.6746 0 3 + 0.99122 6 2 

- 0.63962 9 - 3.4032 (4) 

where 6 = logjo T - 2.1077. These numbers are approx- 
imately 0.7 units of magnitude larger than the corresponding 
values on Earth, which is to say that the same double 



Maximum peak— to— peak = 7.6 microns 



Fig. Dl. Rayleigh wave synthetic seismograms computed for Earth model PREM (bottom) and Mars model AR (top) for similar geometries (A = 90‘ 
azimuth = 20° ) and for equivalent focal depths (30 km), source mechanisms (strike = 45° , dip = 45° , slip 45° ), and moments (M„ = 10 25 dyn cm). 
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(a) 


40 63 


PERIOD (s) (b) PERIOD (s) 

100 158 251 398 40 63 100 158 251 398 




Fig. D2. (a) Logarithmic average excitation of mantle Rayleigh waves for 10 source depths (10-100 km; individual symbols from 0 to 9) for 
Mars model AR. This figure is equivalent to Fig. 3 of Okal and Talandier (1989). (b) Residuals of (a) when equation (4) is used to model the 
excitation. 


couple excites Rayleigh waves with spectral amplitudes at 
the source about about 5 times larger on Mars than on 
Earth. 

The effect of propagation on the spectral amplitude of 
a Rayleigh wave is twofold; g eome trical spreading along 
the sphere, in the amount l/\/sinA, which is independent 
of the radius of the planet, and anelastic attenuation. Once 
again, the latter depends critically on the Q structure of 
the planet, but in any case, the shorter the distances traveled 
generally the less the attenuation. On the other hand, 
Rayleigh wave spectral amplitude is proportional to the 
radius a of the planet (Okal and Talandier, 1989), resulting 
in an additional factor of about 0.5. As a whole, one would 
expect spectral amplitudes to be higher by a factor of a 
few (depending on Q). As shown in Fig. Dl, time-domain 
amplitudes are larger by an even greater factor because 
of significantly less dispersion on Mars. 

Love Waves 

Dispersion characteristics for Love waves are shown in 
Fig. 9. Because of the absence of a low-velocity asthe- 
nosphere in model AR, the Love wave is significantly dis- 
persed, in a fashion somewhat reminiscent of Love waves 
in continental shields on Earth. The excitation of Love 
waves is depicted on Fig. D3. The logarithmic average 
excitability is in general 0.6 units of magnitude higher than 
on Earth and can depend significantly on depth for periods 
less than 65 s. The logarithmic excitation L L can be modeled 
by the best -fitting cubic spline for h = 50 km 


L L = 0.96953 0 3 -0.010307 0 2 

-0.90238 0 -3.1071 (5) 

where 0 = log 10 T - 2.0476. For reasons similar to the case 
of Rayleigh waves, the amplitudes of Love waves are also 
expected to be 10-20 times stronger than on Earth, for 
comparable seismic sources. 

MARTIAN MAGNITUDE SCALES 

In this section we discuss the possible definition of 
magnitude scales for marsquakes. The purpose of using a 
magnitude scale is to relate the size of possible marsquakes 
to our terrestrial experience. Therefore, to convey the sense 
of a source of similar physical size, we must use a magnitude 
scale related to moment, but not to ground motion at 
a given distance, since this concept would itself depend 
on the particular distance chosen. If an earthquake and 
a marsquake shake the ground in the same amount at a 
distance A, the marsquake may be more efficient at 2 A, 
this result being independent of whether A is measured 
in degrees or kilometers. 

Surface Wave Magnitude 

In Okal (1989), we have shown that a theoretical rela- 
tion of the form 


M s = log 10 M 0 - 19.46 


( 6 ) 
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(a) PERIOD (s) 

50 79 128 200 318 



(b) PERIOD (s) 

50 79 126 200 318 



Fig. D3. 


Same as Fig. D2 for Love waves. Equation (5) is used to compute the residuals in (b). 


exists on Earth between the standard 20-s surface wave 
magnitude Ms and seismic moment (in dyn cm), for those 
events small enough that Ms has not saturated. Our pur- 
pose in developing a magnitude scale for Mars will be to 
measure a quantity expected to represent log 10 M 0 - 19.46, 
so that events of similar moments will have similar magni- 
tudes on the two planets. We concentrate in this section 
on the case of the vertical component of Rayleigh waves, 
which are expected to be somewhat dispersed and thus 
may offer the possibility of a time-domain measurement. 

As discussed in the above section, the situation is 
significantly different from Earth because the relationship 
between the apparent time-domain amplitude x(t) and the 
spectral amplitude X(a>) excited by the source is more 
complex than for Earth. In Okal (1989), using phase- 
stationary asymptotics we showed that for 20-s terrestrial 
Rayleigh waves, dispersion was sufficient to “spread out” 
the wave energy over many cycles of the reference period, 
leading to the relationship 


, xwu 

x = 2 

TV^D 


(7) 


where 6k is the range of wave numbers contributing to 
the wave packet, remains significantly larger than 1, which 
in the case of terrestrial Rayleigh waves imposes only the 
loose constraint A > 10°. 

However, in the case of Mars for periods around 20 s, 
the range of k is not significantly different, but the distances 
D are obviously shorter. Most importantly, a is significantly 
reduced (to 0.027 km/s 2 , rather than 0.08 on Earth). 
Consequently, the parameter z reaches 1 only for A > 
50°. 

We will assume, however, that the relationship in equa- 
tion (7) holds, and follow Okal’s (1989) examination of 
the theory underlying the terrestrial Mg. We start with 
the expression 

logio Mq = logio X(a>) + Cg + Cp + 19.38 (9) 

identical to equation (15) of that paper, except for the 
different constant (19.38) reflecting the smaller size of the 
planet. In equation (9) Cg and Cp are source and distance 
corrections, respectively; at T > 45s Cg = -L R , where L R 
is given by equation (4). Phase-stationary asymptotics allow 
the transformation of this expression into 


where D is the distance traveled by the wave along the 
Earth’s surface, U is the group velocity, and a = /dU/dT/ 
is the dispersion parameter. This relation holds as long as 
the quantity 


logio M 0 = log 10 (x / T) + Cg + C D + 0.5 log, 0 A 


+ log,/^^+ 19.38 


( 10 ) 


^5k 


z = 


w 


where 1^ is the length of 1° of arc on the martian surface 
(8) (59.13 km) and a is taken as 0.027 km/s 2 between 17 and 

26 s. At 20 s, the source correction Cg cannot be computed 
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from equation (4), but an independent computation suggests 
the value 2.65 for sources 30-40 km deep. 

The distance correction in the time-domain, Cj D = Cq 
+ 0.5 logio A, could depend significantly on the Q of 20-s 
Rayleigh waves. However, if the martian crust is thick and 
has high Q (in excess of 500), the effect of Q will be 
minimized by the shorter distances traveled, and Cq 0 can 
be modeled by the expression 0.61 logio A - 0.292. These 
values predict that a formula of the type 

Ms = l°gio Y + 0 ’ 61 lo 8io A + 4.21 (11) 

will lead to the correct moment-magnitude relationship 
in equation (6). 

We have tested equation (11) on the synthetic seismo- 
gram shown on Fig. Dl, yielding Mg — 5.77 for a synthetic 
computed with a moment of 10 25 dyn cm. The difference 
of 0.23 units between the computed and expected values 
[according to equation (6)] is probably representative of 
the effect of source mechanism geometry. We also ran simi- 
lar tests at distances of 45° and 135°. At the greater dis- 
tances, the magnitude was found to be very stable (5.81), 
but at the shorter distance, a significantly larger value was 
found (6.56). This illustrates the fact that at the shorter 
distances, dispersion is not sufficient to warrant applying 
equation (7). 

Finally, let us emphasize that the choice of computing 
M s at 20 s was entirely based on an analogy with Earth, 
on which 20-s waves are prominent and well-dispersed. 


When we learn more about the shallow structure of Mars 
and measure actual Rayleigh waveforms, it may be that 
a different period will provide a more justifiable and robust 
measurement. The above formalism could easily be adapted 
to computations at other periods. 

Body Wave Magnitude 

In the case of a short-period body-wave magnitude mj,, 
the proposal of a specific formula relating time-domain dis- 
placement amplitudes to seismic moment is even more 
uncertain. In the terrestrial case, the distance-correction 
term is completely empirical (e.g., Richter, 1958). Between 
distances of 25° and 90°, it oscillates between 6.3 and 
7.0 for shallow earthquakes. Given our rudimentary 
knowledge of the structure of Mars, especially as far as 
attenuation is concerned, we can only surmise that 1-s 
P wave amplitudes can be expected to be approximately 
10 times larger than on Earth at comparable distances and 
propose the formula 

% 

m b-log 10 ^: + 7.6 (12) 

i 1 

(x in micrometers and T close to 1 s). It would then be 
expected that a marsquake and an earthquake of similar 
moments would have have similar body-wave magnitudes. 
Clearly, any further development will have to rely on ac- 
tual teleseismic data from the planet. 
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